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ABSTRACT 



Heat- transfer measurements were made for a single, 
water-heated tube in a pool of R-114 to simulate operating 
conditions of a water chiller. Data were obtained for a 
smooth copper tube, and for two commercially available 
tubes: a spirally roped copper-nickel tube with a porous 
coating; and a copper tube with a structured outer surface 
and a multiple-start helical ridged inner surface. 
Measurements were made for refrigerant-oil mixtures at oil 
concentrations from 0 to 6 mass percent with a boiling pool 
temperature of 13. 8 °C. Results for the two enhanced tubes 
with and without oil are compared to the smooth tube data. 
Enhancement factors for the overall heat- transfer coeffi- 
cient were 4. 0 and 3. 6 for the structured surface and 
porous - coated tubes, respectively, in pure refrigerant and 
at a water Velocity of 2 m/s. For these same conditions the 
enhancement factors for the outside heat- transfer coeffi- 
cient were 14.6 and 6.4 for the porous - coated and structured 
surface tubes, respectively. 
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I. INTRODUCTION 



Recent developments in boiling surfaces have shown 
considerable enhancements in heat transfer performance. The 
worldwide literature on enhanced heat transfer contains over 
3000 published technical papers and reports [Ref. 1: p. 81] . 
This increased interest in heat transfer augmentation is the 
result of incentives for energy and material savings. One 
effective way to improve the heat transfer is by using 
passive augmentation. 

Passive augmentation uses fine-scale alteration of 
surfaces, both external and internal [Ref. 1: p. 82]. These 
surfaces may consist of either an applied porous coating or 
fins deformed in various ways to provide a large number of 
reentrant cavities on external surfaces. Methods of 
internal enhancement include insert devices, forged fins, or 
deformation of the surface (i.e. , spirally roped or corru- 
gated tube). The use of such surfaces can lead to consider- 
able reduction in the size and weight of heat- transfer 
equipment. The reductions may lead to smaller capital costs 
or operating costs or both. 

Of particular interest, in naval applications, is the 
reduction in the size of water chillers in refrigeration 
systems. While a number of investigations are currently in 
progress, the theoretical treatment of the boiling perform- 
ance of various tubes is almost impossible owing to the very 
complicated mechanisms involved. The degree of difficulty 
increases with the presence of oil in the refrigerant 
liquid. In general, refrigeration systems with oil- 
lubricated, hermetically sealed compressors contain small 
mass percents of oil in the refrigerant liquid. Therefore, 
reliable data covering a wide variety of operating condi- 
tions for various refrigerants and different boiling 
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surfaces is essential to further develop more compact and 
efficient evaporators. 

Host of the experiments on enhanced boiling surfaces 
reported in the literature have used electrically heated 
tubes. The difficulties associated with these tubes, espe- 
cially during the instrumentation stage, raise questions as 
to the reliability of the data as described by Wanniarchchi 
et al [Ref. 2: p. 14]. Despite the precautions taken to 
minimize contact resistance, it may still be present between 
the thermocouple locations and the outer boiling surface. 
Another area of doubt, described by Wanniarchchi et al. 
[Ref. 2: p. 14] , is the nonuniform heat fluxes generated by 
commercially available electrical heaters. Additionally, in 
order to use an electrical heated tube, the inside surface 
of the tube must be smooth, or any internal enhancement must 
be bored out. This is disadvantageous as designers would be 
interested in knowing the overall heat - transfer coefficient. 
Most industrial experiments use 2-to 3-meter-long tubes 
with warm water flowing inside them. While the information 
generated from these experiments closely approximate actual 
evaporators, compared to electrical heated tubes, the large 
water temperature drop from inlet to exit (up to 5 K) may 
make it difficult to study the details of the boiling 
process. 

The type of refrigerant selected will also influence the 
design of evaporators. R-114, a moderate-pressure refrige- 
rant, is receiving more attention, in particular for naval 
applications. Advantages to using R-114 include: (a) it 
belongs to the refrigerant group with the lowest toxicity, 
(b) it is very stable with temperature, and (c) it has a 
fairly large value for the energy transfer per unit volume 
of vapor ( i. e. , E v = h (g /v g ) [Ref. 3: pp. 281-291] . The 
advantages of the first two items are clear. The third item 
is advantageous since a higher value for E v means lower 
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pressure drop through the refrigerant piping for a specified 
heat duty. As the normal boiling point of the type of 
refrigerant increases, the E v value increases. R-ll, R-113, 
R-114, and R-22 have E v values in increasing order. R-114 
may be preferable to R-ll and R-113 due to the higher E v 
value and preferable to R-22 due to the ability to use 
lighter components in the refrigeration loop. 

Based on the above discussion, the major objective of 
the present investigation was to study the boiling perform- 
ance of two commercially available tubes, a porous - coated 
tube and a mechanically structured surface, in comparison 
with a smooth tube for the following conditions: (a) tubes 

are water-heated, (b) boiling fluid is R-114 with 0, 1, 2, 

and 6 percent by mass oil, and (c) boiling temperature is 
13.8 °C. 
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II. BACKGROUND 



A. NUCLEATE BOILING 

When the temperature of a boiling surface exceeds the 
saturation temperature of a fluid by a few degrees, nucleate 
boiling occurs. The difference between the boiling surface 
and saturation temperature is the amount of wall superheat. 
Dougherty [Ref. 4: p. 175] defined pool-boiling as vaporiza- 
tion occurring under the following conditions: ( 1) liquid 
depth >> the bubble diameter, (2) there is negligible effect 
on heat transfer due to low or no externally imposed 
velocity of the fluid, and (3) the bubbles move away from 
the boiling surface due to a force field. With the above 
conditions, the heat flow can be expressed as: 

h = -w- ( 2 -d 

A description of the boiling mechanism is given by 
Chongrungreong and Sauer [Ref. 5: p. 701]. A thin layer of 
superheated liquid is formed adjacent to the boiling 
surface. Bubbles nucleate in this thin layer and grow from 
preferred spots on the boiling surface. It is assumed that 
the primary resistance to heat transfer is within the thin 
liquid layer. The height of the liquid in the boiling 
container is not a primary variable, but the liquid proper- 
ties should be controlled. However, large flooded evapora- 
tors may experience a significant "submergence effect" from 
the liquid head. 

3. EXTERNAL SURFACE ENHANCEMENT 

For smooth heating surfaces, bubbles nucleate at various 
scratches and cavities on the surface. Fujii stated that the 
number of active sites increases as the heat flux increases 
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[Ref. 6: p. 48]. Webb [Ref. 7: p. 46] stated that the 
ability to enhance the nucleate boiling coefficient by 
applying some type of roughness has been known for over 50 
years. 

Presently, there are two major types of commercially 
developed roughnesses: (1) porous coating and (2) various 
fins and surfaces with reentrant cavities [Ref. 8: p. 24] . 
The first type of roughness is a sintered particle coating, 
usually copper or aluminum. Figure 2. 1 shows a schematic and 
a photomicrograph of a porous coating. The small reentrant 
cavities are interconnected by substrate tunnels. The 
surface reduces the required superheat for vapor generation 
by entrapping a high density of relatively large vapor 
nuclei in the cavities contained within the porous coating 
[Ref. 8: p. 24]. As described by Webb [Ref. 7: p. 58] , 
researchers found that a critical pore size not a particle 
size governs the number of nucleation sites. Large pores are 
required for liquids with high surface tension and high 
thermal conductivity, while small pores work best for 
liquids with low surface tension and low thermal conduc- 
tivity (i.e. , refrigerants). Czikk and O'Neill [Ref. 9: p. 
53] developed correlations to include the effects of 
different pore sizes of porous- coatings. They concluded that 
there were two resistances to heat transfer for these coat- 
ings: a nucleation superheat related to bubble diameter; and 
a conduction superheat controlled by the liquid film sepa- 
rating the vapor bubbles from the metal matrix. Carnavos 
[Ref. 10: pp. 106-108] found that a porous coating resulted 
in a 700-800 percent better heat - transfer coefficient than a 
smooth surface in a pool of R-ll. Czikk [Ref. 11: p. 98] 
also found increased performance of the porous - coated tubes 
used in his 20-ton water-chiller experiment. 

Structured fins or surfaces form reentrant cavities of 
various geometries. Figure 2.2 is shows an example of a 
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Figure 2.2. Schematic of Manufactured Reentrant Cavity 
( trom Ref. 3 ) . 
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manufactured reentrant cavity. These reentrant cavities act 
as stable nucleation sites, which enhances the heat 
transfer. The nucleation site, to remain active, is depen- 
dent on the mouth diameter falling within a critical range 
and shape with a maximum reentrant angle. This range is a 
function of the fluid properties. The cost of these 
surfaces, presently, are not significantly higher than the 
cost for smooth tubes and they dramatically improve boiling 
performance compared to the smooth tubes. 

The mechanisms by which the reentrant cavities operate 
are complicated. Previous research concluded that the combi- 
nation of bubble evaporation, thermal boundary- layer strip- 
ping and bubble agitation controlled the heat transfer from 
smooth surfaces [Ref. 12: p. 192]. In contrast, different 
experiments conducted by Arshad [Ref. 12: p. 192] and Arai 
[Ref. 13: p. 37] show that thin film evaporation in the 
reentrant cavity controlled the heat-transfer mechanism. 
Bubbles were formed by vapor exiting cavities as the liquid- 
vapor interface of the thin film contacted the cavity 
surface.- Surface tension holds most of the liquid on the 
cavity walls. Ayub [Ref. 14: p. 64] showed this similar 
"thermosiphon mechanism" with his experiments using enhanced 
surfaces . 

Nakayama [Ref. 15: p. 37] gives a fairly detailed 
description of a physical model undergoing this bubble 
growth mechanism. Figure 2.3 shows the three major phases of 
bubble growth. Phase I consists of a pressure buildup in 
the tunnel by evaporation of liquid held in corners and 
continues until the meniscus reaches a hemispherical shape 
of radius r Q = d 0 /2 (where d 0 is the mouth diameter of the 
reentrant cavity). In Phase II, the meniscus grows faster at 
some pores than at others. The Phase I pressure buildup is 
reduced as the vapor enters the growing bubbles. The 
meniscus at inactive pores does not grow due to this vapor 
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Figure 2.3. Physical Model of Bubble Mechanism 
( from Kef. 15 ). 
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pressure reduction. Initially, the bubble expands under the 
high internal pressure, while lacer expansion is controlled 
by the receding liquid inertia around the bubble. Phase III 
is termed the liquid intake phase with liquid flowing into 
the tunnel through inactive pores. This flow is over a short 
interval of pressure depression which occurs as the pressure 
of the bubble and tunnel is lower than the pool pressure. 
The bubble leaves the pore and new meniscus formation closes 
off the pore, ending this phase and returning the cycle to 
Phase I. 

C. EPFECTS OF OIL 

The introduction of oil into the pure refrigerant, in 
general, causes a decrease in the boiling coefficient as the 
oil concentration increases. The effects of oil on the 
mixture are theoretically complicated as most empirical 
equations relating physical properties are derived form 
experimental data and not from theoretical equations. 
Various empirical equations are cited by Chongrungreong and 
Sauer [Ref. 5: pp. 703-705]. Yet, Jensen and Jackman 
[Ref. 16: p. 184] showed that the correlation developed by 
Chongrungreong and Sauer appeared to overpredict as oil 
concentrations increased due to poor prediction of the 
mixture viscosity. 

Predictive equations were developed by Jensen and 
Jackman [Ref. 16: pp. 136-137] for density, viscosity, 
surface tension and specific heat for pure R-113, pure oil, 
and R-113-oil mixtures. These equations are: 



P 



m 




( 2 . 2 ) 




(2.3) 
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(2.4) 



c 




(2.5) 



These equations are based on a physical model shown in 
Figure 2.4 developed by Jensen and Jackman [Ref. 16: pp. 
187-188]. The model shows a bubble growing on a heated 
surface as the refrigerant evaporates from the superheated 
liquid phase into the bubble interior. Due to the less vola- 
tile nature of the oil in the mixture, it does not evaporate 
into the bubble, but, remains at the liquid-vapor interface. 
The decrease in the rate of bubble growth is due to the 
decrease of oil diffusion into the liquid mixture and the 
decrease of refrigerant through the oil layer into the 
bubble. The decrease of the bubble growth rate decreases 
the heat- transfer rate. 

D. INTERNAL ENHANCEMENT 

Internal enhancement increases the heat - trans f er area, 
creates additional turbulence, and produces secondary flows, 
all of which contribute to the increase in the heat 
transfer. The laminar sublayer is assumed to be the prin- 
cipal resistance to heat transfer [Ref. 17: p. 6], The 
developed secondary flow thins the sublayer and moves it 
into the stream where turbulence prevails and without an 
increase in shear; likewise, heat transfer is increased 
without an increase in friction [Ref. 18: p. 30] . Most 
studies show that the friction factor does increase with 
internal enhancement but the amount is strongly dependent on 
the internal enhancement geometry. The increase of both the 
friction and heat transfer is dependent on fin or rib pitch, 
groove depths, and helix angles [Refs. 17,19: pp. 6, 19-21]. 
Additionally, these studies have shown that there are 
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critical geometrical dimensions that give the maximum heat 
transfer. For example, short fins ( << tube diameter) with 
rifling will raise heat- transfer coefficients with no prob- 
lems of flow stagnation. Yet, medium fins may develop flow 
stagnation and decrease heat- transfer coefficients [Ref. 19: 
pp. 22-23] . As such, the selection of internal enhancement 
will generally improve the performance of boiling surfaces. 
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III. EXPERIMENTAL APPARATUS 



A. OVERALL APPARATUS 

A schematic of the experimental apparatus is given in 
Figure 3.1. The photographs in Figure 3.2 show two 
different views of the apparatus. Complete details on the 
design and construction of the refrigerant and oil compo- 
nents are given by Karasbun [Ref. 20: pp. 24-32]. A discus- 
sion of the modifications for the water heating mode is 
provided in this paper. The major components of the appa- 
ratus are: two Pyrex glass tees, an R-114 liquid reservoir, 
a water- ethylene glycol mixture sump, an R-12 refrigeration 
system, a vacuum pump, a water supply tank, three centri- 
fugal pumps, a flow meter, three heaters, a graduated oil 
cylinder, and an oil reservoir. The R-114 boiling and 
condensation occurred in the lower and upper glass tees, 
respectively. The R-114 vapor was condensed by fhe water- 
ethylene glycol mixture pumped through a copper condenser 
coil located in the upper glass tee. A 1/2-Ton R-12 refrig- 
eration system maintained the water- ethylene glycol mixture 
between -18 and -14 °C. 

The major difference between the experimental apparatus 
used by Karasabun and Reilly [Refs. 20,8: pp. 24-32, 30-35] , 
and the one used during the present investigation concerned 
the use of water heating instead of electric heating as a 
means to provide heat to the boiling tube. Filtered tap 
water supplied from a storage tank (15), was pumped by two 
1/2-hp motors (13 and 14) connected in series. A three-way 
valve (V-16) was provided on the downstream side of the 
second pump to obtain two flow paths. The normal flow path 
for data runs was through: the metering valve (V-17), the 
flowmeter (6), the inlet mixing chamber (16), the boiling 
tube, the outlet mixing chamber (16), the chambers provided 
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Figure 3. 1. Schematic of Experimental Apparatus. 
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Figure 3.2. Photographs of Experimental Apparatus: (a) Water 

Inlet Side, and ( L> ) Front View. 





with heaters (7 and 8), and returning to the storage tank 
through valve (V-20). The by-pass flow path was the same up 
to the three-way valve then the path was through heater 
chamber (8) and back to the storage tank through valve 
(V-20). Valves V-17 and V-21 were closed in order to 
isolate the boiling tube. This by-pass option was essential 
for preheating the water to higher temperatures, to achieve 
higher heat flux values in the boiling tube compared to 
values obtained from room- temperature water. A 500-W and two 
1000-W heaters were used to maintain a steady water inlet 
temperature and to preheat the water. Flexible pressure hose 
was used for water piping throughout the apparatus. The 
mixing chambers and sections of hose connecting the mixing 
chambers to the boiling tube were double insulated to ensure 
accurate temperature measurements. 

Copper- constantan thermocouples were used to measure 
thermal emf's for R-114 liquid and vapor, water inlet and 
outlet, and water- ethylene glycol mixture. A series- 
connected thermopile consisting of 10 junctions on either 
end was used to measure the temperature drop of the water 
from the inlet to the outlet of the boiling tube. This ther- 
mopile was calibrated against two quartz thermometer probes 
and agreement was found to be better than ± 0. 02 K. 

B. BOILING TUBE CONSTRUCTION 

Figure 3.3 shows schematics of the smooth tube, and 
two enhanced tubes used. The tubes tested were: 

1. a smooth copper tube; 

2. a 90:10 copper-nickel corrugated tube (commercially 
referred to as Korodense tube) with an external, 
sintered porous coating (i.e. High Flux); 

3. an internally and externally enhanced tube (alloy 
C12200) produced by Wolverine Division of U. 0. P. 

( commerically referred to as Turbo-B); 

4. a porous - coated (High Flux) Korodense tube with the 
porous coating machined off; 

5. a Turbo-B tube with the external enhancement machined 

off. 
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(a) 




(b) 




( c) 



Figure 3.3. 



Schematic of 3oiling Tubes: (a) Smooth Tube, 

( b ) High Flux Coated Tune, and (c) Turbo-3 
Tube. 
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(a) 




(b) 



Figure 3. 4. 



Photographs of Internal Enhancements^ 

(aj Korociense Tube, and ( b ) TurDo-3 Tube. 
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These tubes will be referred to as High Flux, Turbo-B, modi- 
fied High Flux, and modified Turbo-B tubes, respectively, 
with the exception of tubes 2 and 4. Tubes 2 and 4 will be 
referred to as Korodense tubes with respect to discussions 
of only the internal enhancement. All of the tubes had an 
active (i.e. , heated) section of 304.8 mm in length. The 
smooth, High Flux, and modified High Flux tubes had an outer 
diameter of 15. 9 mm and an inside diameter of 12. 7 mm. In 
the case of the High Flux tube, these diameters represent 
the values in the smooth portion of the tube (i.e. , before 
performing the corrugation process). The Turbo-B tube had 
an external structured surface and a multiple- start helical 
ridged inner surface. This tube had an outer diameter of 
16. 9 mm at the base of the external structured surface and a 
minimum diameter of 14. 5 mm at the internal ridge tips. The 
modified Turbo-B tube had identical diameters as the Turbo-B 
tube. The 63.5-mm-long inactive sections on either end was 
insulated by Teflon sleeves of 1. 6 mm wall thickness located 
inside the boiling tubes, as shown in Figure 3.3. The 
requirement for using these two modified tubes is explained 
further in Chapter IV (EXPERIMENTAL PROCEDURES). 

C. COMPUTER- CONTROLLED DATA ACQUISITION AND REDUCTION 

Hewlett-Packard equipment was used for data acquisi- 
tion, control, and analysis, as well as for storing of 
data. An HP-3497A acquisition unit was used to read the 
thermocouple and thermopile outputs. Table 1 gives the 
channel designations used. 

Information entered by keyboard to an HP-9826A computer 
unit prompted and controlled the data acquisition unit. 
Data were analyzed and stored with the same computer unit. A 
step-by-step description of the data- reduction procedure is 
given in Appendix A, along with a printout of the program 
(DRP4) used. 
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TABLE 1 

CHANNEL DESIGNATIONS 



Channel 

0-1 

2 

3 

4 

5 

20 



Function 

R-114 liquid thermal emf's T(l) 
and T(2) 

R-114 vapor thermal emf T(3) 
water- ethylene glycol mixture 
thermal emf T(4) 
water inlet thermal emf T(5) 
water outlet thermal emf T(b) 
thermal emf difference across 
water inlet and outlet T( 7 ) 
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IV. EXPERIMENTAL PROCEDURES 



A. FLOWMETER CALIBRATION 

This section gives a brief description of the flowmeter 
calibration; a more detailed explaination is given in 
Appendix B. A Fischer Porter flowmeter was used to indicate 
percent water flow through the boiling tube. Prior to 
conducting any boiling runs, it was necessary to develop a 
correlation relating the flow percent to the mass flow rate 
and water velocity. Additionally, any temperature effect on 
viscosity would be incorporated into this correlation as a 
correction factor. Calibration runs were conducted over a 
temperature range of 19 °C to 38 °C and varying flowmeter 
settings. The correlation (equation (B.l)) was determined 
from these results. The correction factor for the viscosity 
temperature dependency was determined to be 1. 0 with an 
accuracy of ± 0.02 as explained in Appendix B. In fact, it 
was not possible to find a systematic trend for this correc- 
tion factor with the water inlet temperature. The stated 
accuracy mainly consists of the uncertainties involved with 
the visual reading of the flowmeter settings. 

B. BOILING TUBE AND APPARATUS PREPARATION 

The external surface of the boiling tube was cleaned 
with a Nitol (2% nitric acid and 98% ethyl alcohol) solution 
and both the external and internal surfaces were rinsed with 
acetone and were air dried prior to installation in the 
apparatus. A vacuum test at an absolute pressure of about 50 
mm Hg was conducted on the refrigerant side of the apparatus 
after the boiling tube was installed. If no signs of leaks 
were evident after two hours, the system was pressurized to 
a gage pressure of about 250 mm Hg with R-114 vapor and the 
apparatus was checked for leaks with an automatic Halogen 
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Leak Detector. After fixing any leaks that may have been 
present, R-114 liquid was transferred from the reservoir to 
a pre-determined level in the boiling tee. This level gave a 
free surface 40 mm above the centerline of the boiling tube. 
At the set level, the mass of the R-114 liquid was computed 
to be 2.48 kg. A desired saturation temperature of 4 °C and 
the resultant system pressure were obtained by turning on 
the R-12 refrigeration system and varying the flow of cold 
(approximately -16 °C) water-ethylene glycol mixture through 
the condenser coil with valve V-9. At this saturation 
temperature and pressure, each boiling tube developed 
nucleation sites. The saturation temperature of 4 °C was 
maintained for at least 30 minutes prior to initiating 
nucleate boiling on the complete active length of the 
boiling tube. 

The initiation of nucleate boiling could be done with 
cold water at 19 °C or with warm water at 25 °C through 
valve V-17 or with steam from a steam generation system 
through valve V-22. Early experimental data showed slight 
differences (up to 10 percent) in the computed outside heat- 
transfer coefficient depending on the highest heat flux at 
which the boiling was initiated on each boiling tube. As 
will be discussed in Chapter V (RESULTS AND DISCUSSION), the 
steam intiation method provided the highest initial wall 
superheat and ensured the best possible repeatability of the 
data. This intiation method was used for all of the runs 
unless otherwise specified. 

Steam initiation was conducted by introducing steam 
through valve V-22 and discharging steam drain valve V-18 
for about one minute. After this period, the three-way valve 
V-16 was shifted from a neutral position to the normal flow 
position. Valves V-22 and V-18 were closed and valve V-21 
was opened simultaneously as the three-way valve was shifted 
to the normal position. Valve V-17 was opened wide prior to 
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shifting of the three-way valve. As these valves were shift- 
ed, warm water at 25 ,°C was pumped through the boiling tube 
at the maximum water velocity, which in turn gave the maxi- 
mum heat flux for the given water inlet temperature. 

C. BOILING DATA RUNS 

Once the nucleate boiling initiation was completed, the 
saturation temperature was raised to 13.8 ?C by adjusting 
valve V-9. This valve was also used to maintain this satura- 
tion temperature within ±0. 05 K throughout a boiling run. 
The water inlet temperature was maintained to within ±0. 08 K 
by using the heaters or introducing filtered cold water, 
respectively. The combination of heaters and cold water was 
used to maintain the desired water inlet temperature during 
each run. Two types of boiling runs were conducted at each 
oil mass percent. 

The first type of run was conducted after holding the 
desired saturation temperature, water inlet temperature and 
initial flowmeter setting at steady-state conditions for 
about 10 minutes. After this period, the flowmeter setting 
was decreased with two data sets being taken at each flow 
setting. The two data sets were to show repeatability at 
each setting. The period between pairs of data sets was 
about 5 to 8 minutes, with steady state for each flowmeter 
setting maintained for at least 2 minutes. The second type 
of boiling run consisted of maintaining the previously 
mentioned initial steady-state conditions at a selected high 
water inlet temperature. The run was commenced by intro- 
ducing cold water into the storage tank which lowered the 
water inlet temperature. While holding the water velocity 
constant, six to seven data points were taken, on a one time 
pass, at about every 0. 3 K inlet temperature decrease. When 
the inlet temperature reached a selected low temperature, 
the steady-state conditions of the low inlet temperature, 
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water velocity, and saturation temperature of 13.8 °C were 
maintained for about 10 minutes. Then, the water inlet 
temperature was increased using the heaters, while holding 
the water velocity constant. Again, data points were taken 
on a one time pass, in about 0. 3 K intervals, until the 
inlet temperature reached the previously used high inlet 
temperature. 

Two complete boiling runs of the first type mentioned 
above were performed at 0% and 2% oil concentrations to 
further demonstrate repeatability. Table 2 shows the percent 
difference in wall superheat (i.e. , AT) and heat flux for 
each boiling tube at these two oil percents. The higher 
percentages for the High Flux and Turbo-B tubes for the wall 
superheats are simply because of the low wall superheat 
values of these two tubes when compared to the values for 
the smooth tube. At these low superheat values, a small 
difference in the wall superheat leads to a larger 
percentage difference than that which occurs at higher wall 
superheats. The variation of the water inlet temperature was 
not considered a factor in this percent difference as the 
variation was less than ± 0. 02 K per minute during these 
constant-inlet- temperature , decreasing-water-velocity runs. 

TABLE 2 

BOILING RUN REPEATABILITY 



Tube 


Oil% 


A T 


q 


Smooth 
High Flux 
Turbo-B 
Smooth 
High Flux 
Turbo-B 


0 

0 

0 

2 

2 

2 


4 9 7? 

±30 A 
±30$ 
±10$ 
±10$ 
±30% 


±20% 

±12$ 

Jas! 

r 

± 2% 



Note: initiations performed with steam 

Additionally, data runs were conducted to determine the 
temperature increase across the boiling tube due to the 
pressure drop between the locations of the thermopile 



35 



probes. The data collected showed that this difference was 
less than the ±0. 02 K accuracy of the thermopile and was 
considered negligible. 

Following the boiling runs in pure R-114, oil was 
introduced to the boiling tee through valve V-l from the 
graduated cylinder (5) shown in Figure 3.1. After the 
required volume of oil for a given oil mass percent was 
added, water at the maximum velocity was pumped through the 
boiling tube to promote vigorous boiling. This boiling 
ensured good mixing of the R-114 and oil mixture. Either 
another boiling run was conducted or the system was shut 
down to prepare for another steam initiation. 

D. DATA ACQUISITION AND REDUCTION 

The automatic data acquisition system was prompted to 
record the required thermal emf ' s , by keyboard inputs to the 
computer. The data were immediately processed and printed 
out on a hard-copy printer. A step-by-step procedure of data 
processing and reprocessing is given in Appendix A. Also, 
a printout of the data reduction program, DRP4 , is included 
in this appendix. The heat flux, overall heat-transfer coef- 
ficient, outside boiling heat- transfer coefficient, and wall 
superheat values were computed based on the outside area 
expressed by the diameter if the enhancement of the outer 
surface was removed. Further, in order to account for the 
heat conduction through the inactive end sections, a correc- 
tion was made for heat flux using an iterative computation 
procedure based on natural convection as discussed by 
Karasbun [Ref. 20: pp. 54-56]. 

When a boiling run was completed, the data were repro- 
cessed by computer to obtain the inside and outside heat- 
transfer coefficients for further calculations. A 
Sieder-Tate- type constant ( C j ) for the inside coefficient 
for each of the externally smooth tubes was calculated by 
using a modified Wilson plot [~Ref. 22j . For the smooth tube, 
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modified High Flux, and modified Turbo-B tubes, the outside 
coefficient was calculated simultaneously using the cor- 
relation (equation: (C. 3) and (C. 4)) developed by Rohsenow 
[Ref. 23: p. 969] , and the modified Wilson plot. The 
Rohsenow correlation is based on externally smooth tubes, 
thus the requirement for the modified High Flux and Turbo-B 
tubes. A further explaination of the modified Wilson plot 
is given in Appendix C. The Cj values found for the exter- 
nally smooth tubes were used for the corresponding exter- 
nally enhanced tubes. Knowing the inside coefficient, the 
outside coefficient for the boiling tubes were computed by 
subtracting the inside and wall resistances from the meas- 
ured overall thermal resistance (equation: (C.7)). 
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V. RESULTS AND DISCUSSION 



A. INSIDE HEAT-TRANSFER COEFFICIENT 

Results were obtained for the smooth, High Flux, and 
Turbo-B boiling tubes at a boiling temperature of 13. 8 °C 
and at oil concentrations of 0, 1, 2, and 6 mass percent. 
The heat flux versus wall superheat for these three tubes is 
shown in Figure 5. 1. All three tubes are in the nucleate 
boiling range with decreasing heat flux. It can be seen that 
the High Flux and Turbo-B tubes outperform the smooth tube 
throughout the tested heat flux range. Additionally, the 
Rohsenow correlation for a smooth tube at a boiling tempera- 
ture of 13.8 ,°C with the experimentally determined (from the 
modified Wilson plot) C s f coefficient of 0.0060 is shown in 
this figure. 

The Rohsenow correlation, as used in the modified Wilson 
plot, has an exponent of r = 1/3. Yet, this correlation is 
very sensitive to the r value. A change in this exponent of 
three percent yielded a change in the C S f only of 0.016 
percent but gave a change of 17 percent in the 
Sieder-Tate- type constant inferred from the modified Wilson 
plot. For the purpose of data analysis, the value or r = 1/3 
was used, while knowing that an additional uncertainty was 
being introduced by this exponent. 

For each externally smooth tube, the Sieder-Tate- type 
constant was computed using the modified Wilson plot based 
on 5 to 10 different runs in pure R-114. The computed 
values for the smooth, Korodense, and Turbo-B tubes were 
0.036 ±0.002, 0.066 ± 0.005, and 0.077 ± 0.003, respec- 
tively. For the smooth tube, the experimentally found 
Sieder-Tate- type constant is 33 percent greater known than 
the well-known value of 0. 027 for long internally smooth 
tubes with fully developed flow. This larger constant 
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Figure 5.1. Variation of Heat Flux with Wall Superheat 
at 0% Oil. 



appears Co be the result of the entrance effects of the 
shorter experimental tube and the uncertainty introduced by 
the exponent r of the Rohsenow correlation as discussed 
above. The minimum entrance length for a fully developed 



for Reynolds numbers greater than or equal to 10000. Also, 
as the pipe roughness increases the minimum entrance length 
decreases. 

Withers has developed Stanton number correlations for 
tubes with single-helix and multiple helix ridging 
[Refs. 25,26: pp. 52-56,44]. These correlations are listed 
below. For single-helix ridging, 



pipe flow is given, as a rule of thumb, by Incropera and 
DeWitt [Ref. 24: p. 406] as 



L 



e 
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(5.1) 



D 





(5.3) 



where r and m are determined by tube dimensions, and 



Y 




(5.4) 



For multiple-helix ridging, 





(5.5) 



+ Y 




(5.6) 



r and m again are functions of tube dimensions, 
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A comparison of Stanton numbers was made for the modified 
Korodense and modified Turbo-B tubes using the following 
equations for these two tubes: 




St = 
where 



h. D. 

i i 

Re ?r 




(3.9) 



(5. 10) 



The results for these two modified tubes were compared to a 
copper Korodense (i. e., single-helix ridging) and a 
Turbo-Chill ( i. e., multiple-helix ridging) tubes, respec- 
tively. The latter two tubes tested by Withers were similar 
in radial dimensions but were 1. 5 meter long. The predicted 
St numbers using equations (5.1) -(5.3) for the Withers' 
Korodense and Turbo-Chill tubes were 0.00177 and 0.00273, 
respectively. The experimental St numbers using equations ( 5.9 ) 
and (5.10) for the modified Korodense and modified Turbo-B 
tubes were 0.00242 and 0.00236, respectively. The majority 
of the difference between the predicted and experimental St 
numbers may be attributable to the entrance effects and the 
uncertainity introduced by the Rohsenow correlation. 

3. LIGHT OF? EFFECTS 

As referred to in Chapter IV, Figure 5.2 shows the light 
off effects for the smooth, High Flux (i. e., porous - coated 
Korodense), and Turbo-B tubes. These effects are most prob- 
ably due to the requirement of greater heat flux to initiate 
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Figure 5.2. Light Off Effects at 0% Oil. 



the nucleation sites than the heat flux required to maintain 
the sites once they are activated. As the boiling tube sits 
in the R-114 liquid, prior to initiation, reentrant cavities 
or the nucleation sites flood with liquid. The initial 
formation of the vapor bubbles and thus the initiation of 
the nucleation site requires the input of a greater heat 
flux or higher wall superheat. Figure 5. 3 shows the boiling 
tube initiation with steam for the smooth, High Flux, and 
Turbo-B tubes [Refs. 21,7: p. 478, 62]. As shown in Figure 
5.3, the High Flux tube is more sensitive to light off 
effects. This sensitivity is mostly due to the lower wall 
superheats at which nucleation occurs. A change of . 5 of a 
degree is more pronounced at these lower wall superheats 
than a similar change at the wall superheats for the Turbo-B 
and smooth tubes. 

C. SMOOTH TUBE 

Figure 5. 4 shows the performance of the smooth tube at 
0, 1, 2, and 6 mass percent oil concentration at 13.8 °C. 
Reilly's [Ref. 8: p. 65] data for the electrically heated 
smooth tube at a boiling temperature of -2.2 °C in pure 
R-114 liquid are also included for comparison. The differ- 
ences between the experimental results shown and Reilly's 
data are probably due to physical characteristics introduced 
in the electrically heated tube when soldering the copper 
sleeves (inside of which the electric heaters were fitted) 
on the interior of the tube. As discussed by Reilly [Ref. 8: 
pp. 57-61] , the contact resistance at the interface between 
the sleeve and the inner surface of the tube was minimized 
by tinning. However, the tinning process may not have been 
100 percent successful, thereby introducing unacceptable 
uncertainties into the experimental measurements. Further 
Reilly reported considerable variations (up to 3 K at a heat 
flux of 98 kW/m 2 ) in the measured wall temperatures. Reilly 
attributed this observation to an axially non-uniform heat 
flux generated by the cartridge heater. 
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Figure 5.4. Variation of Heat Flux with Wall Superheat 
for Smooth Tube. 



As shown in Figure 5.4, the wall superheat increased for 
1 percent oil concentration and decreased slightly for 2 and 
6 percent oil. This decrease in the wall superheat for the 
higher oil percentage is possibly due to the enhanced bubble 
formation experienced by the tube because of the foaming 
action of the R-114 and oil mixture decreasing resistance to 
heat transfer as explained by Chongrungreong [Ref. 5: p. 
701] , Nobukatsu [Ref. 27: p. 60] , and Chaddock [Ref. 21: p. 
474] . 

D. HIGH FLUX TUBE 

The performance of the High Flux tube in terms of heat 
flux versus wall superheat is given in Figure 5.5. 
Additionally, this figure shows Reilly's data [Ref. 8: p. 
70] , for the electrically heated High Flux tube in pure 
R-114 at a 6. 7 °C boiling temperature. The presence of oil 
increased the wall superheat, at a heat flux of 40 kW/m 2 , by 
a factor of 1. 22 for 1 percent oil and a factor of 2. 12 for 
6 percent oil. The number of nucleation sites reduced 
considerably and the foaming action increased significantly 
as the oil mass percent increased. These actions are shown 
in Figure 5. 6. With 6 percent oil concentration, it was 
visually observed that the nucleation sites were active only 
in the corrugations, as shown in Figure 5.6 (c). This 
congregation of sites is due, in part, to the sparser porous 
coating on the high points of the tube and the decreased 
wall thickness in the corrugations. As a result of the 
corrugation process during tube manufacture, the wall thick- 
ness at the corrugation is smaller than the rest of the 
tube, which implies lower wall resistance. During the 
sintering process, the coating appears to have concentrated 
in these corrugations giving a thicker coating of copper 
particles and increased number of nucleation sites compared 
to the rest of the tube surface. 
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Fipure 5.5. Variation of Heat Flux with Wall Superheat 
for High Flux Coated Tube. 




U) 




(b) 




The increase of the wall superheat, when oil is present, 
is possibly the result of increased surface tension and 
viscosity. The increasing surface tension and viscosity lead 
to the formation of an oil-rich layer next to the boiling 
surface. This oil-rich layer has an insulating and inhib- 
iting of bubble growth effect within the cavities. The 
higher liquid-vapor surface tension and surface tension in 
the oil-rich layer requires more energy for the formation 
and growth of the bubble. Additionally, diffusion processes 
affect the bubble formation. Oil concentrations near the 
boiling surface are high compared to the concentration in 
the bulk liquid and diffusion occurs from the cavity into 
the bulk liquid. This diffusion continues until an equilib- 
rium point is reached. If oil-rich layers are developing on 
the boiling surfaces, the diffusion from the cavities to the 
bulk liquid slows down aiding the further build-up of the 
oil-rich layer. The vapor bubble forming in the cavity also 
has an oil vapor content which is undergoing its own diffu- 
sion process. As the refrigerant evaporates through the 
oil-rich layer into the interior of the bubble, it must 
overcome large diffusion resistances set up by the oil vapor 
within the bubble. The increased surface tension, viscosity, 
diffusion resistance and the decreasing bubble formation 
rate contribute to the decrease in the heat- transfer coeffi- 
cient and the increase in the wall superheat at a given heat 
flux. [Refs. 5,21,28,29,30: p. 702, 477-478,59,372,82] 

Also, boiling runs were conducted holding the water vel- 
ocity constant and allowing the water inlet temperature to 
decrease and then to increase. The decreasing and increasing 
heat flux with the respective change in the water inlet 
temperature versus wall superheat is shown in Figure 5. 7. 
The motion of the hysteresis ( i. e. , clockwise or counter- 
clockwise) depended on the temperature starting point. If 
the temperature was decreased and then increased, the motion 
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was clockwise. When the temperature was increased and then 
decreased, the motion was counterclockwise. This directional 
change is related to the amount of heat flux required to 
maintain active nucleation sites. In the first case of tem- 
perature decreasing then increasing, less heat flux is need- 
ed to maintain the same amount of sites, but, at the same 
time, heat flux decrease as the water heat capacaity is 
decreasing. In the second case, with the water inlet temper- 
ature increasing then decreasing, the water heat capacity 
increases so the heat flux increases raising the wall 
superheat and the heat- transfer coefficient. 

E. TURBO-B TUBE 

Figure 5. 8 shows a similar relationship as Figure 5. 5 
but for the Turbo-B tube. The increase in wall superheat, 
for a heat flux of 40 kW/m 2 , ranges from a factor of 1.04 
for 1 percent oil to a factor of 1. 37 for 6 percent oil. The 
Turbo-B tube remained active over the full active length of 
the tube, as shown in Figure 5.9. The apparent decrease of 
nucleation sites, as shown by the High Flux tube, did not 
occur with the Turbo-B tube. As discussed pertaining to the 
High Flux tube, the increase in the wall superheat with 
increasing oil concentration is probably due to the 
increased surface tension and viscosity of the mixture 
within the cavities. Also, Boiling runs for constant water 
velocity and varying water temperature were performed on the 
Turbo-B tube and are also shown in Figure 5. 10. The same 
hysteresis motions observed with the High Flux tube are 
observed with the Turbo-B tube. 

F. COMPARISON OF BOILING HEAT-TRANSFER COEFFICIENTS 

The ratios shown in Figure 5.11 are h/h for the 
Turbo-B and High Flux tubes and h/h s for the smooth tube. 
All of the ratios were taken at a heat flux of 40 kW/m 2 . The 
High Flux tube shows the most significant effect due to oil 
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Figure 5.8. Variation of Heat Flux with Wall Superheat 
for Turbo- B Tube. 
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Figure 5. 11. Variation of Boiling Heat-Transfer 
Ratio with Oil Hass Percent for a 
Heat Flux of AO lcW/m 2 . 



concentration which is due to the reasons previously 
discussed. The h/h s ratio for the High Flux tube decreased 
by 65 percent. For the Turbo-B tube the rate of degradation 
of the boiling heat- transfer coefficient is less pronounced. 
The Turbo-B tube shows a decrease of 24 percent at 6 mass 
percent oil. The larger cavity size of the Turbo-B tube, 
compared to the High Flux tube, results in smaller oil 
concentrations in the oil rich layer. This observation may 
explain the smaller degradation experienced by Turbo-B tube. 

G. OVERALL HEAT-TRANSFER CHARACTERISTICS 

The overall heat- transfer coefficient versus decreasing 
water velocity for oil mass percents of 0, 2, and 6 for each 
boiling tube is shown in Figure 5.12. For pure R-114 with a 
water velocity of 2 m/s, the Turbo-B tube outperforms the 
High Flux tube by a factor of 1. 13 and the smooth tube by a 
factor of 4. While the High Flux tube outperformed the 
Turbo-B tube based on the outside heat- transfer coefficient, 
the reverse is true for the overall heat- trans f er coeffi- 
cient. This relationship in U 0 is due in part to the 
increased internal enhancement of the Turbo-B tube. 
Comparison of the Sieder-Tate- type constants for the Turbo-B 
and High Flux tubes (see p. 34) shows that the Turbo-B tube 
has a 17 percent greater inside coefficient than the High 
Flux tube. Both of the enhanced tubes dramatically outper- 
formed the smooth tube. It should be noted that the increase 
in pressure drop due to internal enhancement (compared to a 
smooth- interior case) was not considered here. The pressure 
drop would be an important factor in finally selecting the 
most economical tube type. 
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Figure 5. 12. Variation of Overall Heat -Transfer Coefficient 
with Oil Mass Percent. 



VI. CONCLUSIONS AND RECOMMENDATIONS 



A. CONCLUSIONS 



Considering the data gathered in this investigation for 
boiling of R-114 and oil mixtures at 13.8 °C boiling temper- 
ature, the following conclusions are reached: 



1 . 



2 . 



Based on the overall heat -transfer coefficient, the 
Turbo-B and High Flux tubes outperformed the smooth 
tube by a factor greater than 3, and the Turbo-B tube 
outperformed the High Flux tube by a factor of 1. 13. 



At a heat flux of 40 kW/m 2 with zero mass percent oil, 
the outside coefficients of the High Flux and Turbo-B 
tubes were factors of 14.6 and 6.4 times, respec- 
tively, compared to the smooth tube. These factors 
decreased to 7. 0 and 4. 9 for the High Flux and Turbo-B 
tubes, respectively, with 6 mass percent oil. 



3. The Turbo-B and High Flux tubes showed Sieder-Tate- 
type constants of 0.077 and 0.066, respectively, com- 
pared to the Sieder-Tate constant of 0. 027 for long 
smooth tubes. As noted in the discussion, these values 
may be up to 10 percent lower for internally enhanced 
long tubes than the values for the short tubes. 

4. While the High Flux tube outperformed the Turbo-B 
tube, based on the outside heat- transfer coefficient, 
the High Flux tube is more susceptible to the presence 
of oil. 



B. RECOMMENDATIONS 

Based on the results obtained from this investigation, 
the following recommendations are made: 

1. The investigation of water-heated tubes should be 
upgraded to use long tubes. This would decrease the 
entrance effects that short tubes experience and it 
would come closer to duplicating actual operating 
conditions. Also, the pressure drop effect previously 
mentioned in Chapter V (RESULTS AND DISCUSSION) could 
be investigated in the analysis of the long tubes. 

2. Studies of the boiling heat- transfer coefficients 
should be expanded to tube bundles. The interaction of 
tubes within the tube bundle is a major factor in the 
analysis of heat - 1 rans f er data. The boiling heat- 
transfer coefficient may be strongly influenced by any 
foaming action within the bundle. 

3. Further investigations should be conducted on internal 
tube enhancements due to the major thermal resistance 
being the internal resistance. 

4. Modifications to the condensing side of the present 
apparatus would allow for studies of short water- 
heated tubes at lower boiling temperatures. The major 



58 



modification would involve only the increase 
condenser drain piping to handle increases in 
rate of the refrigerant liquid. 



of the 
the flow 
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APPENDIX A 

DATA REDUCTION PROGRAM 



The data reduction program, DRP4 , used for this investi- 
gation is listed below. A brief description of the major 
sections is given in Table 3, which is followed by an indi- 
vidual line listing of DRP4. 
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TABLE 3 

DRP4 MAJOR SECTIONS 



Line Number 



Description 



10-80 

80-3690 



3700-4165 



4170- 

4225- 

4940- 

6180- 

6415- 

6495- 



4220 

4935 

6175 

6410 

6490 

7280 



7285-7390 

7395-9640 

9645-9725 



9730-9755 

9760-10200 



10205-10360 

10365-10475 



Selection options (i.e. process 
data , plot , etc. ) 

Sub Program Main 
-boiling tube dimensions 
-select desired operating 
conditions (i.e. initiation 
mode, flow rate, Tsat, etc.) 
-convert emf ' s to temperatures 
-compute contact resistance 
(electrically heated mode only) 
-compute various water properties 
-various heat transfer 
calculations( Q , LMTD, U . etc. ) 
-compute various Freon 114 
properties 

-compute natural convective heat- 
transfer coefficient for 
unenhanced ends( iterative 
procedure ) 

-compute heat loss for unenhanced 
ends 

-compute actual heat flux and 
boiling coefficient 
-record and store data 
Functions for curve fits of 
various Freon 114 properties 
Function for polynomial curve fit 
Sub Program Poly 

-calculate nth order polynomial 
Sub Program Plin 
-plot linear graphs 
Sub Program Stats 
-calculate standard deviations 
Sub Program Coeff 
-plot cross-plot files 
Sub Program Wilson 
-calculate modified Wilson 
plot values (see Appx. C) 
Functions for curve fits of 
various water properties 
Sub Program Plot 
-plot non-linear graphs 
Sub Program Symb 
-print various symbols on 
plotted graphs 
Sub Program Purge 
-purge unwanted files 
Sub Program Tdcn 
-calculate water temperature 
rise corrections due to 
pressure drop between 
thermopile locations 
Sub Program U plot 
-print 0 files 
Sub Program Select 
-select various sub 
program options 
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D P r E : Octooer 2 . * ‘.S- 

REVISED : ^ori. j 7. '986 

COM / lap/ lop 
PRINTER IS 1 
CALL Select 

INPUT "WANT TQ SELECT ANOTHER 0P T I0M < 1 = V . 0«N ) ^Msel 
IF I s e 1 = 1 THEN GOTO 40. 

REE D 
BEEP 

PRINTER IS 1 

PRINT "DATA CQLLECTION/REPRQCESSING COMPLETED” 

END 

SUB Tain 
CON /lap/ lap 
COM /Cc/ CC 7 ) , leal 
CON /Mil/ D2 . D i . Do . L . Lu. Kcu 

DIM Emf < 12) ♦ T ( 12) . D 1 a( 9 ) . D2a( 9 ) . D ia( 9 ) , 3oa< 9 ) , La( 9 ) , Lua< 9 ) , Xcua( 9 ) . E t < 19) . 
T nS ( 4 ) [ 1 5 ] 

105 DATA 9. 10036091 ,25727 .34263 . -767345 .3295 . 78025595 .31 
110 DATA -92* 7486589. S.97638E-1 1 ,-2.661 92E-1 3 , 3 . 94073E+1 4 
1 15 READ C(*> 

120 ! DATA "Smoo th" . "H ign F I llx" . "T hernoexe I ~E" . "Thermoexe I -HE" 

1 25 DATA Smooth. Hign FlLLX.TiLrDo-3.Hign Flux Nod.Turoo-3 Nod 

ISO READ Tn5<*> 

125 PRINTER IS 701 
140 BEEP 

f 45 IF Icd-4 THEM PRINTER IS ! 

150 IF Iqp-4 THEM GOTO 1230 

155 INPUT "ENTER MONTH, DATE AND TINE ( MM : DD : HH : NN :$S) " , DateS 
160 OUTPUT 709 : "I'D" : DateS 

165 OUTPUT 709 : "TD" 

-70 ENTER 709: Dates 

175 PRINT 

130 PRINT " Month, date and time Dates 

135 PRINT 

190 P Q INT USING " 1 OX , ""NOTE : Program name : DRP4 

195 SEE 3 

SOU INPUT "ENTER DISK NUMBER". Dn 

205 PRINT USING " 1 6X . ""Disk munoer = "*\ZZ“:Dn 

2 ! 0 BEEP 

215 INPUT "ENTER INPUT MODE ( 0 =3054A . 1 =FILE ) ".In 
220 BEEP 

225 INPUT "SELECT HEATING MODE (0=ELEC; 1 -WATER) Ihm 
220 BEEP 

225 INPUT "ENTER THERMOCOUPLE TYPE ( 0 =NEW . i =OLD > " . Icai 
24Q IF Im=Q THEN 

245 BEEP 

250 INPUT "GIVE A NAME "OR THE RAW DATA FILE" . D2_ r i ieS 

255 PRINT USING "16X.""Ne»j Tile name: "" . 1 4A" : D2_ r i IeS 

260 S A ze I =20 

265 CREATE 3DAT Q2_ r l L eS . S i :e i 

d 7 0 ASSIGN ?F i leL i 0 D2 * i i e ’s 

275 ’ 



30 

35 
4 0 
45 
50 
55 
60 
65 
70 
75 
3 0 
85 
90 
95 
100 



220 1 DUMMY FILE UNTIL Nrun KNOWN 
225 D1 i LeS=”DUMM v " 

290 CREATE 3DAT D t _ f l I eS . S l ze ! 
295 ASSIGN rFile* TO DW’LeS 
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_ 0 0 III J7='JT ;jF • i :Date'5 

305 I” Ihm-0 T HEN 

3 . q 3 

2 T 5 INPUT "ENTER NUMBER OF DEFECTIVE TCS ( O-DEFAULT) ".Id tc 
320 IF Id tc=0 THEN 

325 Late ! -0 

330 Ldtc2=0 

335 PRINT USING ‘’16X,""No defective ^Cs exist 

340 END IF 

345 IF I d t c - l THEN 

350 BEEP 

355 INPUT "ENTER DEFECTIVE TC LOCATION’* ,Ldtc ! 

360 PRINT USING "16X.""TC is defective at location "" . D" : Ldtc 1 

365 Latc2=0 

370 END IF 

375 IF Idtc-2 THEN 

380 SEEP 

385 INPUT "ENTER DEFECTIVE TC LOCATIONS’* . Ldtc !. Ldtc2 

390 PRINT USING "16X.""TC are defective at locations . D . 4X . 0“ ; Ld tc 1 . !_dtc2 

395 END IF 

400 IF Idtc>2 THEN 

405 BEEP 

410 PRINTER IS 1 

415 BEEP 

420 PRINT "INVALID ENTRY- 

425 PRINTER IS 701 

430 GOTO 310 

435 END IF 

4a0 END IF 

445 OUTPUT 2 >Fiiel :Latc! .Ldtc2 
450 ! Im=’ option 
455 ELSE 

460 3EEP 

465 INPUT "GIVE THE NAME OF THE EXISTING DATA FILE" , D2_ r i leS 

470 PRINT USING "16X."*’0id file najne : "" . 1 4A" : D2_f i leS 

4 75 ASSIGN 3»File2 TO D2_fileS 
480 ENTER U>F 1 1 e2 : Nr u.n 

- 1 S 5 ENTER f*F i i e2 : Do 1 dS 

400 PRINT USING "16X."" T his data set taxen on : ! 4A" : Do idS 

495 ENTER *Fi ie2 : Ldtc * , Ldtc2 
500 IF Ldtc!>0 OR Latc2>0 THEN 

505 PRINT USING " 1 6X , ’’"The mocoup les were defective at locat i ons : "" . 2< 3D . 4X ) 

Latcl .Ldtc2 

5 1 0 END IF 

515 ENTER $File2:Itt 

520 END IF 

525 ! IF Ira-0 AND Ihm= ! THEN T 595 
530 3EE D 

535 INPUT "VifiNT TO CREATE A PLOT FILE 7 ( Q-M . ! =Y >’Mp lo t 

540 IF Ip lo t= 1 THEN 

545 BEE? 

550 INPUT "GIVE NAME FOR PLOT FILE" . P_* 1 ie$ 

555 CREATE 3DA 7 P__ r iie3.4 

560 ASSIGN 2 *P Lot TO P_f iie'5 

565 END IF 

570 IF Ihm= I THEM 

575 5EE D 

530 INPUT "WANT TO CREATE Jo FILE 7 ( 0 = N . ; - O ’’ . Inf 
535 IF I af = 1 THEN 

530 BEEP 
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o 0 0 

■SOS 
S ' 0 
6 15 
820 
625 
630 
835 
6*0 
645 
650 
655 
660 
665 
670 
675 
680 



INPUI 

CREA! 



:gn 

END IE 
BEEP 
INPUT “' 
IF Ire- 
BEEP 
INPUT 



- :j * - - 
3DA “ 

i'ij f : : 



JO 



■ i . j r u*3 



IAN 7 IQ 
THEN 



I.REATE Re FILE'? ( 0-N . ! -v ) , ; 



ENTER Re FILE 



NAME 

CREATE 3DAT Regies. 10 
ASSIGN i'Refiie TO RefiieS 
END IF 
END IF 
PRINTER IS 1 
IF In-0 THEM 
BEE° 

PRINT USING “4X 
IF Ihm-0 THEN 



.Ref ile5 



Select tube number' 



635 


PRINT 


USING 


“5X . 


““0 


smooth 4 inch Ref 




690 


PRINT 


USING 


“ 6 X . 




Smooth 4 men Cu (Press 


/Slide) 


695 


PRINT 


USING 


“SX . 


.. t.o 


Soft bolder 4 men uu 




700 


PRINT 


USING 


“ 6 X . 


..MO 


Soft Solder 4 inch HIGH 


FLUX 


705 


PRINT 


USING 


“SX , 


““4 


'Aieiand Hard 3 inch 




710 


PRINT 


USING 


“SX . 


““5 


HIGH FLUX 3 men 




715 


PRINT 


USING 


“SX . 


““6 


GEUA-< IS Fins/in 




720 


PRINT 


USING 


“SX . 


....? 


GcHA-K 2S Fins/m 




725 


PRINT 


USING 


“SX , 


m“3 


GEwA- T 19 r i ns/ i n 




730 


PRINT 


USING 


“SX. 


'• «« Q 


GcHA- T 25 F ins/ in 




725 


ELSE 












740 


PRINT 


USING 


“SX . 




Smooth tube 




745 


PRINT 


USING 


“SX. 


..... 


H i q n r 1 ux 




750 


PRINT 


USING 


“SX. 


•M.O 


Turoo-3 




755 


PRINT 


USING 


“SX. 


.» ».o 

, O 


High Flux Hod 




760 


PRINT 


USING 


“SX, 


•••• 4 


Turoo-B Mod 





765 
770 
775 
730 
735 
730 
755 
30 0 
305 
3 I G 
315 
320 

gnc 

330 

335 

840 

345 

350 

355 

360 

365 



““Tube Number: 
“ T ube NumDer : 
“Tube Type: 



" . D“ : I tt 
, DD“ : 1 1 1 
. I 5 A” ; TnS< 1 1 1 ) 



END IF 
INPUT Ttt 
OUTPUT Fiiel : It: 

END IF 

PRINTER IS 701 

IF 1 1 1 < 1 0 THEN PRINT USING “ 1 6X , 

IF 1 1 1 > 9 THEN PRINT USING ”1 6X . ' 

IF Ihm* I 'HEN PRINT USING ”1 6*X , ’ 

3E~ D 

INPUT “ENTER OUTPUT VERSION < O-LONG . 1 -SHORT. 2-NONE) “. lov 
BEEP 

INPUT “SELECT ( 0=L IQ . ! WAP . 2= < LIQWAP > /2 > “.I lev 

DIMENSIONS FQ "ESTED TUBES 
ELECTRIC HEATED ^ODE 

D 1 ®D l ame rer at thermocouple positions 

DATA .011 1 125 . .IT I 1 125 . .01 11 125 , .0123540 . .012446. ,0129540 . .01 0096: 
DATA .0100965. .01157. .01 157 
READ D1a(*> 

Q 1 =D 1 a< 1 1 1 > 



870 ! D2-Diameter of test section to the base of fins 

375 DATA .015375. . 0 1 5375 .. 0 1 5375 . .015324. .015875. .015824, .01270 

330 DATA . 0127. .0133. ,0138 

335 READ D2a< - > 

390 ! 

395 ! Di-Instde diameter of unennancea ends 



64 



905 READ J i a< - » 

9 1 0 ! 

915 ! Do-Outside diameter of unennanced ends 

920 DATA .015375, .015375. .015875. .015324. .015375. .01 5324, .01270 . .01 270 01331 . 

.01321 

925 READ Doa(-) 

930 ! 

925 I L=Lengtn of enhanced surface 

940 DATA .1016. .1016. .1016, .1016. .2032, .2032. .2032. .2032. .2032, .2032 

945 READ La(-) 

950 ! 

955 ! Lu=Length of Linennanced surface at the ends 

960 DATA .0254. .0254. .0254. .0254. .0762. .0762. .0762. .0762, .0762. .0762 

965 READ Lua(-) 

970 ! 

975 ! <cu-Thernal Conductivity of tube 

980 DATA 398.344.344.45.344,45,344.944.398,398 

985 READ Kcua<^> 

990 IF Ihm-1 THEN 

905 » 

1000! Data statements for ?jater heating mode 
1005! 

1010 DATA 0.015875.0 .015875,0.0169.0.0138.0 .0169.0.0 .0.0.0 
1015 READ D2a( * > 

1 020 DATA 0 .0127.0.0127 .0.0145,0.0127,0.01 45.0 .0 .0.0 ,0 
1025 READ D ia< * > 

1030 DATA 0.015875.0.015375,0,0169,0.015875.0.0169.0.0.0.0 .0 
1035 READ Doa< - > 

1040 DATA 0.3048.0.3043.0.3048.0.3048.0.3048,0 .0.0.0 .0 
1045 READ La<-> 

1050 DATA 0.0254.0.0254,0.0254.0.0254,0 .0254,0 ,U. 0.0.0 
1 055 READ Lua( ♦ > 

1060 DATA 398 . 45 . 398 . 45 . 398 . 0 . 0 . 0 , 0 . 0 
1065 READ Kcua<*) 

1070 END IF 
1 075 D2 = D2a< 1 1 1 > 

' 080 D 1 s D ia( I tt > 

1 085 Do=Doa( 1 1 1 > 

1090 L='_a( I tt > 

1 095 Lu=Lua( 1 1 1 > 

1 1 00 Kcu=<cua( 1 1 1 > 

1 1 05 Xn a . 3 
1110 p r=.3 

1115 IF It t-0 THEM Cf-l.70E-9 
1120 IF 1 1 1 > 0 THEN Cf -3 . 7037E+1 0 . 

1 125 A=PI*< Do *2-Di 2)/4 
1130 P= p I -Do 
1135 BEE P 

M4Q INPUT -tube INITIATION MODE. ( 1 “HOT VIA TER .2-STEAM. 3-COLD HATER)’* . It im 

1145 IF Itim =1 THEN PRINT USING '* 1 6X . r u£>e Initiate: Hot Water * 

1150 IF Itim*2 THEN PRINT USING " 1 6X . ‘’’' x uDe Initiate: Steam”"*' 

1155 IF Itim-3 THEN PRINT USING ’* 1 6X . ’*** nice Initiate: Cold Water'**’" 

[’60 INPUT “TEMP/VEL MODE: ( 0= T-CGNST . V-DEC :!- T-QEC . 7-CONST : 2-T-INC , V -CONST) ” . 
I tv 

1165 IF Itv = Q THEN PRINT USING *’ 1 6 X , *’’*Temp/Vei Mode: Cons tan t/Dec reas 1 ng‘ 

1170 IF I tv - 1 THEN PRINT USING "1 6 X . ”*'Temp/Ue 1 Mode: Decreas 1 na/Constant ””” 

'175 IF I tv = 2 THEM D RINT USING "1 6 X . ""Temp/Vei "*ooe: Increasing/Constant 

1 1 30 INPUT ‘"WANT TO RUN WILSON PLOT'? ( 1 = Y . 0-N > “ , Iw 1 i 

1 1 35 IF Ihm= T AND I ijj 1 1=0 THEN 
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\ - ^ ~ k w i_ 

“95 IF lit- 1 OR It— I 'HE:: Ci-.OSr 
1 200 IF Itt-2 UR Itt-* rLi EiJ Ci* . 062 
12G5 3EE P 

T 2 1 G INPUT ’‘ENTER Cl < DEF : NH- . 032 . HF= . 059 . TB= . 062 )". C: 

1215 PRINT USING " 1 SX . ""3 ieaer- 'are 

1220 PRINT USING "16X."" Constant = . Z . *D ,# :Ci 

1225 END IF 

1220 IF Ihm«» AND In- 1 AND luiii*! THEN 
1235 IF 1 1 1 = 0 THEM Ci=.032 

1 240 IF Itt-f OR Itt-3 THEN Ci=.U59 

1245 IF Itt-2 OR Itt*4 THEN Ci=.062 

1250 ASSIGN ?F-le2 TO - 
1255 CALL W 1 1 son ( Cf . C i > 

1250 ASSIGN *Fiie2 "0 D2_ r iie5 

1255 ENTER i le2 : Mr an . DoidS .Ld tel , Ldtc2 . I t t 

1270 END IF 

1275 Nslld = 0 

1230 IF Iop- 4 THEN Ihn=! 

1235 IF Ihm-1 THEN Nsud-8 

12S0 N tc*b 

1295 J=i 

1300 Sx-0 

1305 Sy = 0 

1310 Sxs=0 

1315 S x y = 0 

1 320 Repeat : ! 

1325 IF Im-0 THEM 
1330 Dt ld-2 . 22 
1335 Iao*2 

1340 ON KEY 0.15 RECOVER 1320 
1345 PRINTER IS 1 

1350 PRINT USING "4X . ’“‘SELECT OPTION 

1355 PRINT USING "SX . ""0-TAKE DATA 

1360 IF I h m = 0 THEN PRINT USING "SX . = SE" HEAT -LUX""" 

1365 IF Ihm-1 THEN PRINT USING "SX.""1=SE7 HATER FLOW RATE 

1370 PRINT USING "SX.""2=SE7 Tsat 

1375 PRINT USING "4X . ""NOTE : KEY 0 = ESCAPE 

1380 BEEP 

1385 INPUT Ido 

1390 IF Iao>2 THEN Ido =2 

1 395 IF Ido-0 THEN 2' 45 

1 400 ! 

14Q5! LOOP TO SET HEAT FLUX OR FLOWMETER SETTING 

1410 IF Ido-1 THEN 

1415 IF Ihm-0 THEN 

1420 OUTPUT 709 ; " AR AF’2 AL 1 3 VR5" 

1425 BEE D 

1430 INPUT "ENTER DESIRED Odo".Daao 

1435 PRINT USING "4X . ""DESIRED Gao ACTUAL Qdn 

1441] E r r - 1 0 0 0 

1445 FQR >1 TO 2 

1450 OUTPUT 709; "AS 3A" 

1455 S un = 0 
1460 r OR J t - T TO 5 
1485 ENTER 709 :E 
1470 S um = S urn * t 
^ 4 75 NEXT Ji 

1480 IF I-» r HEM Vo it -Sum-5 
1485 IF 1=2 THEN RnD=t 
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* -so * i e :< t : 

’ 435 ^qgp*^o : **^no / ( ° I ~~J2 * L > 

1500 IF ABS< Aqap-Qoap > >Er r THEN 

1505 IF Aacp>Daap THEN 

’510 8EE D 4000, .2 

1515 3EEP 4000. .2 

1520 BEEP 4000. .2 

1525 ELSE 

1 530 SEE* 250 . .2 

! 535 SEEP 250.. 2 

1540 BEEP 250, .2 

1545 END IF 

1550 PRINT USING ”4X . MZ . 3DE .2X . MZ . 2DE” : Qacp , AqdP 

1555 WAI T 2 

1560 GOTO 1445 

1585 ELSE 

1570 BEEP 

1575 PRINT USING ”4X , MZ . 3DE , 2X , MZ . 2DE” : Dadp . AqdP 

1530 Er r = 50 0 

1535 HA IT 2 

1590 GOTO 1445 

1595 END IF 

1600 ELSE 

1605 3EEP 

1610 INPUT '‘ENTER FLOWMETER SETTING”. Fms 
1615 GOTO 1350 
1620 END IF 
1625 END IF 
1620 ? 

1635! LOOP TO SET Tsat 
164Q IF Ido = 2 THEN 
’645 IF Ikdt = 1 THEN 1670 
1650 3EE? 

1655 INPUT '‘ENTER DESIRED Tsat”. Dt id 

1660! PRINT USING “4X.”” QTsat ATsat Rate Tv Rate 

1665 Ikdt=I 
1 670 Old 1 = 0 
1675 01a2=0 

1630 Nn= 1 
1635 Mrs=Nn MOD 15 
1630 Nn=Nn+i 
1695 IF N rs = 1 THEN 

1700 PRINT USING ”4X.”“ Tsat Tidt T ld2 Tv Tsuitid TinLet 

Tout 

1705 END IF 

1710 IF INm-0 THEN OUTPUT 709:“AR AF3 ALU VR5“ 

1715 IF Ihm= i THEN OUTPUT 709:”AR AF0 AL5 UR5” 

1720 FOR 1=1 TO 6 

1725 IF lhm®0 AND I>4 THEN 1360 

1730 S i_Lm = 0 

! 725 OUTPUT 709: "AS SA” 

1 740 FOR Jl-i tq 20 
1745 ENTER 709:Eiiq 
1750 Sujn = Sum+t 1 io 
1755 NEXT Ji 
1 760 El iq = Sum/20 
1 765 Tld^NTvsvi El ia) 

'770 IF 1=1 THEN r ld!*Tlc 

’775 IF 1=2 THEN Tla2-ia 

•730 IF 1 = 3 THEN ~v-id 



Tp i le 
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T35 IF I-J 'sur)d = " la 

*790 i F I - “D HEN ~ t n t e r = ' i c 

i 795 IF I “ *j HEM ^out-^id 

1800 NEXT I 

1805 IF Ihm= • THEM 

1810 OUTPUT 70S : ”AR AF20 AL20 VR5" 

1315 OUTPUT 709: "AS SA" 

1320 Sam = 0 

1325 FOR <'< = ' '0 20 

1330 ENTER 709 :E 

1835 S urn = S tun -t 

1840 NEXT <Ic 

1345 Emf < 7 > =ABS < Sam/20 > 

1850 Toile-Emf <7)/3.96E-4 
1855 END IF 

1360 Atld«< Tldl +Tid2>*.5 

1885 IF AB$( Atia-Dt ld> > ,2 THEM 

1370 IF A t id>D t Id THEN 

1375 BEEP 4000 . .2 

1330 BEEP 4000, .2 

1885 BEE 3 4000, .2 

18S0 ELSE 

1895 BEEP 250. .2 

1900 BEEP 250. ,2 

1905 BEE 3 250. .2 

1910 END IF 

1915 Err 1 -At id-0 Id 1 

1320 0 id 1 -At id 

1925 Err2-Tv-Qla2 

1930 0id2-Tv 

1935 IF lhm = 0 THEM 3 R INT USING ”4X . 5 < MOD . DD . 2X ) M : D t id . T id 1 . T ld2 . Tv . Tsump 
1940^ IF Ihm= 1 AND ldp = 0 THEN PRINT USING ”4X . 7 ( MDD . DD , 2X ) " ; D t id . T Id 1 . T Ld2 , Tv . Ts 
loud . T t n iet . T d t ie 

■945- IF _Ihm= 1 AND Idp-4 THEM PRINT USING "4X . 5 (MDD . 3D . 2X> . 3 ( M3D . DD . 2X) M ; D t id , T 1 
dl .TId2.Tv.7gLunp.Ttn 
7 950 HAI T 2 
1955 GOTO 1635 
I960 ELSE 

1965 IF ABS ( A t id-D t id) > . ! THEN 

1970 IF Atid>Dtid THEN 

1975 BEE 3 3000 . .2 

1930 BEEP 3000 . .2 

1985 ELSE 

1990 BEEP 900 , .2 

1995 3EEP 300. .2 

2000 END IF 

2005 c.rri=Atid“'3Ia1 

201 0 0 ia 1 =A t id 

2015 Err2=Tv-0ia2 

2020 0 Id2 = T v 

2025 IF Ihm-0 THEM PRIMT USING ”4X . 5< NDD . DD . 2X ) " : 0 t Id . Tid 1 . T La2 . Tv . Tsumo 
2030- IF Ihm= 1 ^ THEN PRIMT USING "4X . 5< HDD . DD . 2X ) . 3< M3D . DD . * X ) M : D t ic , T Id 1 . T Ld2 . Tv 
, ”samD , r iniet.Tpiie, 



2035 


HA I ^ 




2040 


GOTO 


r 635 


2045 


ELSE 




2050 


BEEP 




2055 


Er r 1 


= A t* 1 d - 0 1 d 


2060 


0 la I 


id 


2065 


Err2 


= T v-i][d2 
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2070 Oia2- 7 v 

20 75 : r Ihm-0 ~~FN PRIN 7 JSING ’4* . 5< k lDD .20 . 2X > " : D 1 1 d . ~ Id i . T ld2 . 7* . Tsotid 

2080 IF Ihm-1 THEN PRINT USING ”4X .SCMDD . DO .2X>" : Dt ic . T Id 1 . TId2 . T v . Tsunc . 7 ime 

, Id i le . Tout 

2085 WAIT 2 

2090 GOTO 1685 

2095 END IF 

2 t00 END IF 

2105 END IF 

2110! ERROR TRAP FOR Ido OUT OF BOUNDS 

2115 IF Ido>2 THEN 

2120 3EEP 

2125 GOTO 1250 

2130 END IF 

21351 

2140! TAKE DATA IF In-0 LOOP 

2145 IF Ikoi=1 THEN 2165 

2 1 50 3EEP 

2155 INPUT "ENTER BULK OIL r\3oo 
2160 Ilcoi-1 

2165 IF Ihm-Q THEN OUTPUT 709:"AR AFO ALM VR5 M 
2170 IF Ihm-1 THEN OUTPUT 709;'*AR AFO AL5 VR5” 

2175 IF Ihm-Q THEN Ntc«f2 
2180 FOR 1-1 TO Ntc 
2185 OUTPUT 709 : "AS SA" 

2190 Sum-0 
2195 FOR Ji-1 TO 20 
2200 ENTER 709:E 
2205 Sutm-$um+E 

2210 IF I - < 9-Msud> OR I=<10-Nsud> THEN Et<Ji-l>-E 
22 1 5 NEXT Ji 

2220 Kd 1=0 

2225 IF 1 = < 9-Msud > OR I-dO-Nsub) THEN 
2220 Eave-Sum/20 
2235 Sun-0. 

2240 FOR Jk *0 TO 19 

2245 IF ABS< E t ( Jfc ) -Eave) <5 . OE-6 THEN 

2250 Sum-Sum*E t ( Jk ) 

2255 E ! SE 
2260 Kdl-<di*1 

2265 END IF 
2270 NEXT Jk 

2275 IF I - ( 9-Nsud ) OR I-(!0-Nsud> THEN PRINT USING "4X.""Kdl - . DD" : Kdl 

2230 IF Kdi> 1 0 THEN 
2235 BE^p 
2290 BEE 0 

2295 PRINT USING "4X.""Too fnucn scarier ina in data - repeat aata ser’‘ ,,,, 

2300 GOTO 1345 

2305 END IF 
2310 END IF 

2315 Enf(I)-Sum/<20-Kdl> 

2320 NEXT I 

2325 IF Ihm-f T HEN 

2330 OUTPUT 7Q9:"AR AF20 AL20 VR5 ” 

2335 OUTPUT 709 : "AS SA" 

2340 S._Lm = 0 
2345 FOR :<k = i TO 20 
2350 ENTER 709 :E 
2355 Sum-Sum*;! 

2360 NEXT Ki^ 

2265 inf ( 7 ) = ABS < Sum > / 20 
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2375 IF Inn = ‘j ~HEN 

2330 OUTPUT 709 : **^R AF : 2 AL ' 3 VR5" 

2335 PGR I 3 ' TG 2 
2230 OUTPUT 709 : "AS SA” 

2335 Sum-0 
2400 FOR Ji-I ~G 2 
2405 ENTER 703:2 
2410 Sum^Sum^t 
2415 NEXT Ji 

2420 IF I-t THEM Vr-Sum/2 
2425 IF 1 = 2 THEM Ir = S.jun/2 
2*30 NEXT I 
2435 END IF 
2440 ELSE 

2445 IF lhm = 0 • THEN ENTER W i Ie2 : Sod , T o 1 dS . Emf < * ) . Ur , I r 
2450 IF Ihm=1 THEN ENTER *F i ie2 : Bop , To ids . Emf ( * > , Fms 
2455 END IF 
2460! 

2465! CONVERT emf'S TO TEMP . VOLT . CURRENT 

2470 Twa=0 

2475 FOR I=i TO Mtc 

2480 IF I d t c > 0 THEN 

2485 IF I-Ldtc! OR I-Ldtc2 THEM 

2490 T ( I >*-99.99 

2495 GOTO 2545 

2500 END IF 

2505 END IF 

2510 IF 1 1 1 < 4 AND THEN 

2515 IF I > 4 AND I<9 THEN 

2520 T<D— 99.99 

2525 GOTO 2545 

2530 END IF 

2535 END IF 

2540 T ( I ) -r NTvsv < Ent < I ) > 

2545 NEXT I 

2550 IF 1 1 1 <4 THEM 

2555 FOR I=* TG 4 

2560 IF I=Ld tc 1 OR I=Ldtc2 THEM 

2565 Tua=Tija 

2570 ELSE _ 

2575 T ua=T)ja i *T ( I ) 

2580 END IF 
2535 NEXT I 
2590 Tija=Tua/ ( 4- Id tc ) 

2535 ELSE 

2600 IF Ihm= 1 THEM 2650 

2605 FOR I=i tq p 

2610 IF I = L d t c ! OR I=Ld tc2 THEN 

2515 T ua= T ya 

2620 ELSE _ 

2625 T ua= T wa* T < I > 

2530 END IF 
2625 NEXT I 
2S4Q T lja = Tua/ < 8- la tc ) 

2645 END IF 
2650 TId=T(S-Msuo) 

2655 Tld2 -T( iq-Msud) 

2660 TIda= ( Tld^Tid ^; - . 5 
2665 Tv = T ( * 1 -Nsud ) 
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2675 


~ ia2=-39 . 99 




2630 


r v«< 7 < : n ) + 7 < ; ; 


) >/2 


2635 


END IF 




2630 


IF lhm=0 THEM 


2710 


2695 


Tsijjnp = T < 1 2 -Mslld } 


2700 


Tin let* I ( ! 3-Ms 


i in ) 


2705 


T out=‘(" 4 -Msud 


> 


2710 


IF I h m = 0 THEN 




2715 


Amp=A8S < I r ) 




2720 


Volt*ABS<Vr>*2 


5 


2725 


Q=Uo 1 t^Amp 




2730 


END IF 




2735 


IF Itt-0 AND I 


hm = 0 


2740 


Kcu=FNKcu< Tua) 




2745 


ELSE 




2750 


Kcn s <cua( 1 1 1 ) 




2755 


END IF 





2760! 

2765! FOURIER CONDUCTION EQUATION NITH CONTACT RESISTANCE NEGLECTED 
2770 IF Ihm-0 THEN Tu-Tua-Q-LOG< D2/D1 ) / < 2-P Wcu*L ) 

2775 IE Ilav-0 THEN Tsat=Tlaa 

2730 IF Ilqv-1 THEN Tsat- < T lda-Tv ) * . 5 

2735 IF I i civ =2 THEN Tsat-Tv 

2730 IF Ihm* 1 THEN 

2735 7avQ= T !niei 

2300 Graa~37 . 9353+ . 1 Q4388~Tavg 

2305 TdroD*ABS< Emf < 7 ) ) ♦ \ . E + 6/ ( 1 0»Gr ad) 

2310 Tavgc^Tin ier-Tdrop* . 5 

2315 IF ABS( Tavg-Tavgc) > . 0 1 THEN 

2820 Tavq= ( Tavg^Tavgc ) ♦ . 5 

2325 GOTO 2300 

2830 END IF 

2335! 

2340! COMPUTE HATER PROPERTIES 
2345 Kw-FNKu< ^avg) 

2350 Muua= H MMuu < T avg ) 

2855 Cpw^NCdu < T avg ) 

2360 Pru*»’NPr , j< Tavq) 

2865 Rhow =f: NRhoijj< Tavg ) 

2370 Tui-Tavg 
2875! 

2380! Compute MDOT 

2335 Mdot*3.9657E-3+Fms-<3.S1955E-3-Fms-<3.320Q6E-S-Fms*< I . 23633E-7-Fms-4 . 3 1 397 

E-i 0 ) > ) 

2390! Mdo t-Mdo t*< 1 . 0365- 7i n 1 e t*< 1 . 96644E-2- 7 in let*5 . 252E-6 > ) / 1 .0037 
2395 Kat-0 

2300 Q*Mao t*Cp'd*Tarop 

2905 Lmtd=Tdrop/LQG( < T iniet-Tsat)/ < T miet-Tdrop-Tsat) > 

2910 Uo=Q/ < PI *Do*L < *Lm td ) 

2915 Ru = Qo~L0G< Do/Di ) / < 2 . *<ca) 

2920 Tu=Tsat td 

2925 Vw-Mdo r/ < Rhow*^ 1^0 i 2/4) 

2930! IF Kdt=0 ^HEN 
2925! Kdt-i 

294Q! "drop® Tdrop- . 004^V W • ^ 

2945! GOTO 3710 
2950! END IF 

2955 Reu®Rhou-Wu*0 i /Muua 

2960 H i =C * D i *Reu . ru ( 1 /3 . ) ■* * Mu.ua/ FMMuu < Tu i ) ) .14 
2965 njic-7avg-0/<PI*Qo-L-^i ) 
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29 70 IF 4BS< J ! ~nEN 

.1975 ; = < Tij i *!■ r u i c > * . 5 

esso .sot a :?bo 

2935 END IF 

2390 Tui-<Twi^Tu iC )^.5 

2995 Ho- 1/<1 /Uo-Do/<0i*Hi >-Ru> 

3000 ThetaD=Q/ ( Ho* p I*Do*'_ ) 

3005 Tu-Tsat-Mhetao 
3010 Thetao- Tu-7sat 
3015 IF Thetab<0 THEN 
3020 BEEP 

3025 INPUT "TUALLCTSAT (O-CQNTINUE. 1 -END )'Mev 

3030 IF Iev-0 THEN GOTO <325 

3035 IF lev- 1 THEN 3535 

3040 END IF 

3045 END IF 

3050! 

3055 ! COMPUTE VARIOUS PROPERTIES 
3060 Tf 1 1m 3 < Tu+Tsat ) * . 5 
3065 Rho 3( - MRho< T f i Im ) 

3070 fi iL- FNMu( Tf 1 1 rn ) 

3075 K-^-MK C Tf i Im> 

3030 Cp-FNCp<Tf ilm) 

3085 Be ta 3> - NBe ta< Tf i lm) 

3090 Hf q=FNHf a ( Tsat) 

3095 Ni=Nu/Rho 
3100 A 1 pna=</ < Rho^Cp ) 

3105 Pr-N i/A ipna 
3110 Psat 3 FNPsat ( Tsat > 

3115! 

3120! COMPUTE NA T URAL-CCNVECTIVE HEAT- TRANSFER COEFFICIENT 
3125! FOR UNENHANCED END($> 

3130 Hbar* I 90 

3135 Fe= < Hbar^P/ ( Kca^A ) ) ' . 5-*Lu 
3140 Tanh=FN Tanh < Fe ) 

3145 Thera 3 The taD*-Tann/ Fe 

3150 Xx“ < 9 . 3 t *Be ta* The tao^Do 3^ T ann/ < Fe^N l ♦A lpna) ) .166667 
3155 Yy = ( 1 * < . 55'3/Pr ) (9/16)) (3/27) 

3160 Hbarc=*VQo-< . 6- . 337 • <x/ Yy ) 2 

3165 IF ABS< < Hbar-Hbarc) /Hbarc > > . 00 1 THEM 

3170 Hbar = < Hbar +Hbarc > ♦ . 5 

3175 GOTO 3135 

3130 END IF 

3185! 

3130’ COMPUTE HEAT LOSS RATE THROUGH UNENHANCED END(S) 

3195 Q 1 = ( Hbar *P*<cil*A ) . 5 * The t ao - i anh 
3200 Qc-0-2*OI 

3205 As 3 PI~D2-L 
2210! 

32*5! COMPUTE ACTUAL HEAT FLUX AND BOILING COEFFICIENT 

3220 Qap-Oc/As 

3225 H tiLDe=OdP/ The tao 

3230 Cs f - < Cp • The rao/ Hf ) / ( Qdp/ < Mu-Hf n ) - < . 0 1 4/ ( 9 . 3 1 -^ho > .5> 

?'^c I 

3240’ RECORD TIME OF DATA TAKING 
3245 IF lm=0 THEM 
3250 OUTPUT 709:" T D‘* 

3255 ENTER 709:ToldS 
3260 bND IF 
3265 » 



M/3. )*Pr 
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::70' r o winter 

3275 »RINTE* 13 "0: 

3230 IF Io«/-0 ~HEM 
7235 Q 9I f !7 

3290 PRINT USING "T OX . ""Data Set Nuuimer - . DDD .2X . '"'3uik Gi. % = "" . DD . D . 5X . 1 

4ft'' : J . 3 od . o Id'S 
3295 IF lr,m=0 THEM 

3300 PRINT USING " 1 OX . '• ,, TC No: 12 3*557 

1 

3205 PRINT USING ”1 OX . "“Temp : ,,M . 3 ( : X . MDD . DD ) " : T< 1> , T < 2 > . T (3 ) . T< 4 ) , T ( 5 ) . T < 6 ) , T ( 
7 > . T < 3 > 

3310 PRINT USING "lOX."" T ja Tliqd ^iiQd2 Tvapr Psat Tsump 

3315 PRINT USING M 1 OX . 2(MDD . DD . 1 X) . IX .MDD .DD . 1 X .2( 1 X .MDD . DD) . 2X , MDD . D" : Tua.TId. 
T ld2 . Tv , Psat . Tsump 

3320 PRINT USING "lOX."" Thetas Htube Qdo 

3325 PRINT USING ”1 OX .MDD .3D . 1 X .MZ . 3DE . 1 X .MZ . 3DE” : Thetao .Htooe . Q op 
3330 ELSE 

3335 PRINT USING "lOX."" Fms Vu T S at Tinl Tdrop ThetaD q Uo 

Ho 

3340*_PRINT USING M ! OX ,4 < 2D . 3D . 1 X > . Z . 3D . : X . DD . DD . 1 X . 3 < MZ . 3DE . I X > '* ; Fms , Vw . Tsar . T i 
nlet. Tdrop. Thetao. Gd 
3345 END IF 
3350 END IF 
3355 IF Iov-1 THEN 
3360 IF J»1 THEN 
3365 PRINT 
3370 IF Ihm-0 THEM 
3375 PRINT USING “ 1 OX . RUN No 
3330 ELSE 

3385 PRINT USING ,, 10X. ,,M FMS 
3330 END IF 
3395 END IF 
3400 IF Ihm-0 THEN 

3405 PRINT USING " T 2X . 3D . 4X . DD . 2X . MDD . DD . 3 ( 1 X . MZ . 3DE) ,f ; J . Boo . Tsar . H tuoe . Qdp . The 
tao 

3 4 1 0 “ L U E 

3415 PRINT USING " • 2X . 3D . 4X . DD . 2X .HDD . DD . 3( 1 X . MZ . 3DE ) " : Fms . 3 op . Tsat . Htooe . QdP . T 
he tao 
3420 
1425 
3430 
3435 
344 0 
3445 
3450 
3455 
3460 
3465 
3470 
?4?5 
3480 
3485 
3490 
3495 
3500 
3505 
35 : 0 
3515 
3520 
3525 
3530 



END IF 
END IF 

IF In-0 THEM 
BEEP 

INPUT "OK TO STORE THIS DATA SET < 1 -Y . 0-N) ?” .Ok 
END IF 

IF Ok - 1 OR Ira-1 THEN J-JM 
IF Ok-i AND Im-0 'HEN 

IF Ihm-0 THEN OUTPUT 3F 1 le i ; B op . T oiaS . Emf ( ♦ ) . Ur . I r 
IF Ihm-1 THEN OUTPUT '*F l ie 1 ; Bop . To ids . Emf < ♦ > . Fms 
END IF 

IF Inf- ! THEN OUTPUT 5IJf 1 ie ; U.j . Uo 
IF I re- 1 THEN OUTPUT ^Ref Lie: r ns.Reij 

IF ( Ira- 1 OR Qk - ! ) AND Ip lot -I THEN OUTPUT PP 1 or : Oao . The tab 

IF lm»0 r HEM 

BEEP 

INPUT "MILL THERE BE ANOTHER RUN < 1 =/ . 0-N ) ?". Go_on 
Nr an- J 

IF 0o_on=0 THEN 3535 
lF Go_on< > 0 THEN Repeat 
ELSE 

IF J<Nrun^ f THEN Repeat 
END I- 



Gii/l 


Tsat 


Htooe 


Qdp 


1 he tap 


OIL X 


TSAT 


HTUBE 


QDP 


THE^AB'” 
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7.535 IF I*-} THEM 
35-0 3EE~ 

3545 -'PINT US I MG " ’ OX . “"NOTE : ““.Z2. ,, ' , lara runs 'jere scored m r 1 le ““ . ! QA“ : p- 
. D2_ ; : ieS 

3550 ASSIGN i>F iia\ TG - 
3555 OUTPUT i»F 1 ie2 :Nrun- * 

3560 ASSIGN ?<Filei TO D T f iie‘5 

3565 ENTER i le 1 : Dare‘5 . tc 1 . LatcZ . 1 1 r 
3570 OUTPUT PF 1 ie2 : DateS . Ldtc i . Ldtc2 , 1 tt 
3575 FOR 1-1 TO Nr.un-1 
3530 IF Ihm-0 THEM 

3585 ENTER ? J F i le I : Sod . To IdS ,Ent( * > . Vr . Ir 
3590 OUTPUT 3Fi ie2 : Bop . ToidS .Emt < * > . Vr . Ir 
3535 ELSE 

3600 ENTER ^F 1 le I : Boo . To IdS . Emf < *> . Fms 
3605 OUTPUT *F 1 le2 ; 3 od , T o IdS . Emf <*>, Fms 
3610 END IF 

3615 NEXT I 

3620 ASSIGN *F 1 le i TO * 

3625 PURGE “DUMMY” 

3630 END IF 
3635 BEEP 
3640 PRINT 

3645* IF Iplot-1 THEM PRINT USING “ * 0X . ““NOTE : ““.ZZ.““ X-Y oairs jere stored m 
Dior data file ““ . 1 
3650 ASSIGN >Fiie2 '0 - 
3655 ASSIGN *Piot TO * 

3660 IF Inf-T THEM ASSIGN ?Ufile TO - 

3665 IF Ire-1 THEN ASSIGN *Refiie TO * 

3670 CALL Stats 
3675 BEEP 

3630 INPUT “LIKE TO PLOT DATA ( 1 -Y . 0-N ) T\ GV 
3685 IF QV = 1 THEM CALL Plot 
3630 SUBEND 

3695 ' 

3700! CURVE FITS OF PROPERTY FUNCTIONS 
3705 DEF -NKcad) 

37 1 0 ! GFHC COPPER 250 T 0 300 < 

37 1 5 TV = T + 273 . 15 !C TO K 

3720 K- 434 -. ll2*rk 

3725 RETURN K 

3730 FMEMD 

3735 DEF FNMa(T) 

3740? 170 '0 360 < CURVE C IT OF VISCOUSITY 

3745 TV = T+273 ♦ 15 !C TO K 

3750 Mu® EXP < -4 . 4636^< 1 0 M . 47 /TV ) ) * 1 . 0E-3 

3755 RETURN Mu 

3760 FMEMD 

3765 DEF FNCo<T) 

3770! 130 r 0 400 < CURVE FIT OF Cp 
3775 TV = f +9 73 . 15 • 0 T P < 

3730 Cp= . 40 1 33-*’ l* . S5007E - 3* T V +1 .5 I494E-6- TV 2-6 . 67353E- ! 0* Tk 3 

3785 Cp=Cp-U)00 

3790 RETURN Cp 

3795 FNEMD 

3300 DEF FNRho< T ) 

3305 TV = T + 273 . !5 'C ^0 ** 

3310 :<=!-< 1 . 3* f V /7 53 . 95 ) f K M3 R 

3315 Ro®-6. 32-61 . 14b4t4-< u V3> ► « 6 . 4 i 30 1 5* < - » 7 . 4 76333* < .5+? . : • 9323* X '2 

J320 R o = p o / . 062423 
3825 PE' URN Ro 
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9830 FNEND 
28?5 DE r "NPr<T> 

33*0 Pr — NCs< T ) -- NMll( T > /F NK ( T ) 

3845 RETURN D r 
3350 FNEND 
3355 DEF FMK ( T ) 

3360! T <360 < NITH T IN C 
3365 K = . 071 - . 000261 *T 
3370 RETURN .< 

3875 FNEND 

3380 DEF FNranh(X) 

3385 P*EXP< X) 

3890 Q-1/P * 

3395 Tann«<P-0)/(P+0> 

3900 RETURN Tann 

3905 FNEND 

3910 DEF FNTvsv(V) 

3915 COM /Cc/ C(7) .leal 
3920 T=C(0> 

3925 FOR 1=1 TO 7 
3930 T-T+C(I)^ri 
3935 NEXT I 
3940 IF leal- I THEN 

39*5 T«T-S.7422934E-2*T*<9.02770*3E-3-T*< -9 . 32599 1 7E-5 ) ) 

3950 ELSE 

3955 T-T+8.62G897E-2 + TW3.76199E-3-T-5.0639259E-5) 

3960 END IF 
3965 RETURN T 
3970 FNEND 
3375 DEF FNBeta(T) 

3930 RoD=-NRho< T+ . 1 > 

3985 Rom=i-NRho( T- . t ) 

3990 Be ta=-2/ ( Rop+Rom ) ♦< Rod- Rom ) / . 2 

3995 RETURN Beta 

4000 FNEND 

4005 DEF FNHf g< T) 

4010 Hf g= i .3741344£>5-T-<3.309436IE-2 >T-i .2165143) 

4015 RETURN Hfg 

4020 FNEND 

4025 DEF FNPsat(Tc) 

4030! 0 TO 30 aeo F CURVE FIT OF Psat 
4035 Tf = 1 . 3 * t c +32 

40 4 0 P a= 5. 345525 -**Tf ♦< . i 53520 82+ Tf ♦< i ..434Q96oE-3 + Tf*9.bl50671E“6) ) 
4045 Pg= p a- 1 4 . 7 

4050 IF Rg>0 THEN ! * = p $IG.-=m Hg 

4055 Psat = p g 
4060 ELSE 

4065 3 -ar= p g-29.32/l4.7 

4070 END IF 
4075 RE'URM p sat 
4030 r NEND 

4035 DEF FNHsmoo tn < X . 3 op , I sat > 

4090 DIM 4(5) .3(5) ,C<5> .D<5> 

4095 DATA .20526. .25322. .319048. .55322. .79909. 1 .00253 
4100 DATA .74515. .72992. .73139. .7*225. .68472. .841*7 
4105 DATA .41032. . ’7726. .25142. .54306. .31916. ! .0845 
4110 DATA .71403. .72913. .72565. .696691 . .665367. .61839 
4' 15 READ A ( * ) , 3 ( * > . C < ♦ > , D < * > 

4120 IF 3op<6 THEN I = Bop 
4125 IF 3op=6 THEN 1=4 
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-r :o 
41 35 

4 ! 4|J 
41 4S 

4 150 
4155 
4180 
4165 



6495 
6500 
6505 
6510! 
6515 
5520 
6525 
6530 
5535 
6540 
6545 
6550 
655! 
6555 
6560 
6555 
6570 
6575 
6580 
5581 
6535 
5590 
6595 
6600 
6605 
6610 
6615 
5620 
5625 
6630 
6635 
564 0 
6645 
6650 
6655 
6660 
6665 
657 0 
6675 
6680 
6685 
o630 
o635 
6700 
6705 ' 
6710 
67 1 5 
6720 
6725 
6730 



IF 3od = 0 1 = 5 

IF Isa:=’ ~HEM 
Hs=EXP< A< I)+B< D-lGGCO ) 
ELSE 

Hs-tXP(C< I)*0< D-LQGCO ) 
END IF 
RETURN 6s 
FNEND 



SUB 4 1 Ison < Cf . Ci ) 

COM /Wii/ D2 , D i . Do . L . Lil . Kcu. 

DIM Emf < 12) 

4LIS0N PLOT SUBROUTINE DETERMINE CF AND CI 
BEEP 

INPUT "ENTER DATA FILE NAME" .FiieS 
BEEP 

PRINTER IS 1 



PRINT USING **4X 



Select ootionL 



PRINT 


USING " 4X . "** 


0 Vary 


Cf 


ana 


C i 


PRINT 


USING "4X." M 


1 Fix 


Cf 


Vary 


Cl 


PRINT 


USING "4X . "" 


2 Vary 


Cf 


Fix 


Cl 


PRINT 


USING "4X . "** 


3 Fix 


Cf 


Fix 


Cl 


INPUT 


"ENTER OPTION* 


* . Icf 1 X 









PRINTER IS 


701 










IF 


Icf L X = 0 


THEN 


6535 








IF 


Icf ix>0 


THEN 


BEEP 








IF 


I C f 1 X = 1 


THEN 


INPUT 


"ENTER 


Cf 


. Cs f 


IF 


Icf ix = 2 


THEN 


INPUT 


’•ENTER 


r r 

s, 1 


. C i 


IF 


Icf i x = 3 


THEN 


INPUT 


"ENTER 


Cf . 


Cl" 



. Csf 



PRINTER IS 



INPUT "Nan t To Varv Coer f 7 ( ’ = Y . 0=N ) Mccoet 
IF Iccoet-1 THEN INPUT "ENTER C0EFF".R 
PRINTER IS 701 

IF Icfix-Q OR Icfix-2 THEN Cfa-.004 

IF Icf ix = ' THEN Cfa-Csf 

Cia-Ci 



Xn- .3 



rr a ,j 
Jj-0 
R r= 3 . 

IF Iccoef = 1 THEN Rr = R 
PRINTER IS * 

PRINT Do.Di „<ca 
ASSIGN *F:ie TO FiieS 

ENTER <F i le : Mr u.n . DateS . Ldtc ! . Ld tc2 . 1 1 1 
Rw = Do*l 3G ( Do / D i > / ( 2 ** cu.) 

3 x = 0 
Sy = 0 
S x y = 0 
$ x 2 = 0 
Sv2=0 

FOR 1=1 TO Nr an 

ENTER i'F i ie : Bod . To iaS . Enf <• > . Fns 
CONCERT EMF '3 TO "EMPERATURE 
FOR J-! T 0 5 
T( J ) =*- NTvsv ( Emf ( J ) ) 

NEXT J 

Tsar = < T ( 1 ) + T ( 2 ) >*.5 
Tavg = T ( 5 > 
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3 / 5 j l' 3.C = !< / » 9 8 5 O ^ . ■ i 4 O O ■* 1 av Q 
674Q Tcrop-Enf ( 7 ) * ' . E t ^/ ( T 0 . ♦uraa > 

6745 7 avcc= T ( 5 >-7arop** . 5 

6750 IF ABS< Tavg-Tavgc ) > . 01 i HEN 

6755 7avq= < Tavgt*Tavgc> ♦ . 5 

6760 GOTO 6735 

6765 END IF 

6770! 

5775! Compute properties of water 
6730 Kw-FNKw ( Tavg ) 

5735 Muua = FNMuu ( Tavg > 

6730 Cpw^FNCpw ( T avg ) 

5795 Pru*FNPrw< Tavg > 

6800 Rhou = FNRhouU Tavg ) 

6305! 

6310! Compute properties of Freon-114 

6315 Lmtd-Tdrop/L0G<<T(5)-T sat )/<T<5)-Tdrop-Tsat> > 

6320 IF J j = 0 THEN 
6325 Tw-Tsat+P’r^um td 
6330 Thetab-Tu-Tsat 
6335 Jj=1 
6340 END IF 
6845 Tf -< Tu+Tsat ) * .5 
6350 Rho-FNRho<Tf > 

6355 Mu=FNMu< T f ) 

6860 K-FNK(Tf) 

6865 Cp=FNCp< Tf ) 

6370 Beta-P’NBe ta( Tf ) 

6875 Hf g=FNHf a ( Tsat ) 

6380 Ni-Mu/Rho 
6835 A lpna=</ ( Rho*Co > 

63S0 Pr-Ni/Aipha 
6895! 

5900! Analysis 
6905! COMPUTE MOOT 
6910 A= P I*< Do 2-9 i '2)/4 
6915 P= p I*Do 

6920 Mdot-3 . 9657E-3+t* ns*< 3 .6 1 955E-3-Fns^< 8 . 32006E-6- p ms*< 1 .23638E - 7- p os^ ii . 3 1 397 
E- 1 0 > > ) 

6925 Q = Mdo t^Cp<j*TaroP 

6930! COMPUTE NATURAL -CONVECTIVE HEAT- TRANSFER COEFFICIENT 
6935! FOR UNENHANCED END<3> 

694Q Hbar^iSO 

8945 Fe- < Hbar ♦P/ < Kcu^A ) > ' . 5-*Lu 
6950 T an n = FH T an n ( Fe > 

6955 Theta=ThetaD*Tann/Fe 

6960 Xx= < 9 . 3 1 ♦Be ta~Tbe tao^Do 3*7anh/ (Fe*N i ♦A l^na) > . 1 66667 

6965 Yy = < 1 *< .559/Pr) " <9/1 6) > (8/271 

5970 Hbarc=</Do*< . b 1 * . 337^<x/ Yy > '2 

6975 IF AB$ <( Hbar -Hbarc > /Hbar )>. 0 0 1 THEM 

6930 Hbar = < Hbar^Mbarc ) ♦ . 5 

6985 GOTO 6945 

6990 END IF 

6995 ! 

7000! COMPUTE HEAT LOSS RATE THROUGH UNENHANCED ENDS 

7005 01 * ( Hbar* p *-\cu~A ) . 5-TnetaD-Tann 

7010 Oc-0-2-01 

70 1 5 As= p I~D2-L 

7020! COMPUTE ACTUAL HEAT FLUX 

7025 Qcid = 0c/As 
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703 0 
7035 

7 040 

7045 

7050 

7055 

7060 

7065 

7070 

7075! 

7030 

7085 

7090 

7095 

7100 

7105 

7110 

7115 

7120 

7125 

7130 

7135 

7M0 

7 1 45 

7150 

7 1 55 

7160 

7165 

7170 

7175 

7130 

7185 

7190 

7191 

7 1 95 
7200 
7205 
72" 0 

72 • 5 
7220 
7225 
7220 
7235 
7240 
7245 
7250 
7255 
7260 
7255 
7270 
7275 
7230 
7235 
7290 
7295 
7300 
7305 
7310 
7315 
7 320 
7 3 7 5 



IF lc*.x*0 IR I c f . a = 2 HE; J .sr = \Cf < * . / R r > 

The t ao=Cs f /Cp^hT q- < Qdp / ( Hu.-*’* * q ) * < . 0 " 4/ ( 9 . 3 1 ^Rho > > .5) 1 1 /R r ) r 
Ho = 'jap/ The taD 
0meqa=Ho/Cf 
U o = 0 / < P I ♦Oo^L *Ln td ) 

Vy = Mdo t / < Rhoy* p I*-D i '2/4) 

Rey ®Rhoy ~Uy-»-D i /Ht±ya 
Tuj i = Tu^Q~Pw/ (PI *Do^L ) 

Gama-Ky/D i -*Pey . 3* p ry < 1 /3 . ) *< Muya/FNMuy ( Ty i ) ) '.14 

PRINTER IS 1 

Yy=< 1 . /Uo-Ry ) ♦Omega 

Xy-Qmeqa*Do/ ( Gama~D i > 

Sx^Sx-’-Xw 
Sy = Sy Yy 
Sxy-$xy + Yy*Xy 
Sx2~$x2+Xy *-Xu 



Sy2=Sy2+ Yy~ fy 
NEXT I 

ASSIGN »F lie TO - 

M- < Sx*$y-Nrun*Sxy ) / < Sx*Sx -Nrun-*Sx2 > 

C = < Sy-Sx^M ) /Nrun 

IF Icfix-0 OR Ictix-3 OR Icfix-4 THEM 

Cic-t/M 

Cfc=1/C 

END ^ IF 

IF let i x = 1 THEM 
Cic-I/N 
Cfc-Cf 
END IF 

IF Icfix-2 THEM 



C jlc=G i 



Cfc-1/C 
END IF 

IF Icfix=3 THEN 7280 

IF ABS< <Ci-Cic>/Cic>> .001 OR AB$< < CF-Cf c ) /Cfc) > .001 THEN 
Ci 3 < C ; + C ic > - . 5 
Cf*<C:+Cfc)*.5 
PRINTER IS 1 

PRINT USING ”2X . *** Csf = . NZ . 3DE . 2X . '*** Ci - ’*'* .MZ .3DE” ;Csf .Ci 

PRINTER IS 701 
GOTO 6655 



END IF 
PRINT 

PRINTER IS 701 

PRINT USING ’*23X . **'* Cf Ci 

PRINT USING "3X . "“ASSUMED . MZ . 3DE . 3X . MZ . 3DE" : C f a. C : a 

POINT USING "3X . ‘"“CALCULATED , \M . 2DE . 2X . MZ . 3DE” :Csf . Ci 
PRINT 

Sun2*Sy2-2^~Sxy-2*C*Sy 4 -M"2*Sx2*2^ f 1*G-Sx- l “Mrun*0 '2 
PRINT USING ”1 OX . ”*'S»un or Squares = , Z . 3DE" : Sluti2 

PRINT USING -1 OX. -•Coefficient = ,D . 3DD" ;Rr 



SUBEND 

DEF r NNu.y < T ) 

A-247.3/< Ti-133. 15> 

Nu=2.4E-5-M] 'A 
RETURN ‘"a 
FMEMD 

DEF FNCoy < T ) 

Cpw- 4 .21 120853- 2. 2S326E- 3- T^< a. *42361 E-5^2. 71 423E-7-T) > 

RETURN Coy* 1000 
FNEND 
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7330 
1 ^>o5 
7340 
7345 
7350 
7355 
7380 
7365 
7370 
7375 
7380 
7385 
7390 



10265 
1 0370 
1 0375 
1 0380 
1 0385 
10390 
10335 
1 0400 
1 0405 
10410 
10415 
10430 
10425 
10430 
10435 
1 0440 
10445 
1 0450 
1 0455 
10460 
! 0465 
1 0470 
10475 



3EF - "iRhotj < 7 > 

Ro = 999. 52946-- < . 0 1 269- >< 5 . 4825 1 2E-3- T * I .224147E-5) > 
RETURN Ro 
F MEND 

DEF FNPru ( T> 

P ru 3 *- NCpu < i ) *»■ NMuw < T ) /FMKu < T > 

RETURN Pru 
FNEMD 

DEF FNKw < T ) 

X-<7+273. T5)/273 . 15 

Kw = - . 92247+X*< 2 . 3395 -X*< 1 .300 7 -Xm . 52577- . 07344-X ) ) ) 

RETURN Kw 

FNEMD 



SUB Select 
CON /IdD/ Idp 
BEEP 



PRINTER IS • 
PRINT USING 
PRINT USING 
PRINT USING 
PRINT USING 
PRINT USING 
PRINT USING 
PRINT USING 
PRINT USING 
PRINT USING 
INPUT Icp 



M 4X. '^Select option: 

”6X.””0 Taxing data or re-processing previous data 
, *SX." M T Plotting data on Log-Log 

"6X. ,,,, 2 Plotting data on Linear 

"6X . "**3 Make cross-plot coer'ft file'*"" 

"6X. M ”4 Re-c i rculate water 

,# SX.”*'5 Purge 

”SX . '*’*6 ^-Dtod correction 

•’SX.-Y Print Uo File 



IF 


Icd = 0 


THEN 


CALL 


Main 


IF 


Ido= 1 


THEN 


CALL 


Plot 


IF 


IdD = 2 


THEM 


CALL 


Pirn 


IF 


lap = 3 


THEN 


CALL 


Coef 


IF 


IdD = 4 


THEN 


CALL 


Mai n 


IF 


Iap=5 


THEN 


ChLl 


°urg 


IF 


Iap = 6 


THEN 


CALL 


T den 


IF 


Iap = 7 


THEN 


CALL 


Uopr t 



SUBEND 
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APPENDIX B 

FLOWMETER CALIBRATION 



A Fischer Porter Flowmeter was used in the experimental 
apparatus to indicate the water flowrate. Prior to 
conducting the boiling tube data runs, calibration of the 
flowmeter was performed over a temperature range from 19 °C 
to 38 °C. The goals of the calibration were: 

1) to develop a mathematical expression relating flowm- 
eter percent reading to mass flowrate and 

2) to determine a viscosity correction factor due to 
varying temperature range and flowrates. 

A weigh tank, platform- type scale, and stopwatch were used 

to record water collection data. 

Data were processed using the program "FMCAL" , listed 
later in this appendix. This program took the flowmeter 
percent reading, water weight collected, and elapsed water 
collection time as the inputs. The outputs were an experi- 
mental mass flowrate and a mass flowrate difference. This 
difference was given by a mass flowrate computed from a 
correlation minus the experimental mass flow rate for the 
same flowmeter percent reading. The correlation was based on 
a fourth-order least- squares fit to the calibration data. 
After all the data runs were completed, the data points were 
combined into one file and a single fourth-order polynomial 
was generated to describe the relationship between the 
flowmeter percent reading and the mass flowrate. Of the 
original data points, 98.61 percent were within a ±5.5 
percent range. Reprocessing the data points led to 78 
percent within a ± 1.02 percent range. The equation (B. 1), 
resulting from the calibration data analysis, is a fourth- 
order polynomial from a curve fit of the reprocessed data 
without a viscosity correction factor: 
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m 



(3. 1) 



= 5.96 37e- 3 - 3 . 1690e - 3 ( FMS) 

-3 . 8 0 2 0 e - 6 ( FMS ) ~ * 1 . 2369e- 7 f FMS) 
-4 . 3l90e- 10 (FMS) 4 

(FMS = flowmeter percent reading) 



Experimental discharge rates are less than the rates 
calculated from theoretical flow equations primarily because 
of the effect of internal molecular friction of the fluid 
and the viscosity. The theoretical flow equation is: 



n 




(3. 2) 



where Q=volumetric flowrate, A =cross-sectional area of the 

w 

narrowest part of annulus, V,=volume of float, Pc= density 
of float, A-=cross-sectional area of the largest part of 

float, and o..=density of fluid. The Coefficient of 

w 

Discharge, C . , is a function of viscosity. A constant value 

G 

for in equation (3.2) would be desirable as this would 

G. 

show a negation in the variations due to the viscosity. The 
effect of the float shape on C . with varying Reynolds number 

• Cl 

is shown in Figure 3.1 (a)’. The Reynolds number for the 
flowmeter uses an effective diameter of the tube inner diam- 
eter minus the float outer diameter. This plumb bob style 
float presents a large surface area to the fluid stream and 
the viscous effects increase as the flow rate increases. 
[Ref. 31: ?p. 9304-9808] 

The data collected from the calibration runs showed 
little or no change in the Discharge Coefficient over the 20 
to 90 percent flowmeter setting range and the 15 K change 
in fluid temperature. The square edge plumb bob float 
reaches a constant Coefficient of Discharge at Reynolds 
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(a) 




(b) 



Figure 3. 1. Effect of Float Shape on Coefficient 
Discnarce: (a) Snuare Edce Plumb Bob 

^trom Ref. 30), and (b) Taper Edge PI 
r 1 o a t . 
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of 

Float 

umb 



number > 10000. This is due to the effective shape of the 
float on the streamline pattern as shown in Figures B. 2 (a) 
and (b). The experimental calibration range of 1300 < Re < 
27000 displayed a flattening out of the Coefficient of 
Discharge. A different effective shape is seen by the 
streamline pattern because of the minimal clearances between 
the tube bead guides and the float. Also, the float used in 
this particular flowmeter had tapered edges versus square 
edges. This different effective shape is shown Figures B. 2 
(c) and (d). Figure B. 1 (b) shows the effect of the float 
shape on the Coefficient of Discharge for the tapered edge 
float . 

Additionally, the accuracy of the flowmeter contributes 
to the flattening effect on the Coefficient of Discharge 
shown by the experimental data. The accuracy effect is felt 
in two major components, reproducibility and scale factor. 
Reproducibility accuracy is minimally determined by the 
change in the fluid flowrate, corresponding float movement, 
and observation of float position. The calibration data had 
good reproducibility over the various flowrate and tempera- 
ture ranges. The scale factor accuracy is inherent to the 
scale markings on the flow tube. Fischer Porter conducted 
numerous tests of flowmeter calibration readings. They used 
a maximum scale fraction error of 0. 9 mm on a 250 mm tube. 
Based on their results, Fischer Porter reported [Ref. 31: 
pp. 9814-9816] an accuracy of plus or minus one percent at 
high flow rates was a reasonable expectation. Also, at low 
flow rates ( < 20 percent) a resonable expectation would be 
an accuracy of eight to ten perecent. 

The experimental data followed a trend of negligible 
viscosity effect. The flow tube design, float shape, and 
scale factor accuracy were the main contributors to this 
negation of a viscosity correction factor. As such, equation 
(B.l) was used without a viscosity correction factor to 
calculate the mass flowrate for the boiling tube runs. 
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(b) 



Figure Q. 2. Effect of Float Shape on Streamline Patt 
( a) Square Edge Plumb Bob Float ( trom Re 
and ( b ) Taper Edge Plumb 3ob Float. 



84 



H.fD 



APPENDIX £ 
MODIFIED WILSON PLOT 



As referenced in Chapters III and IV, this appendix 
provides the data reduction procedures used to compute 
inside and outside heat - trans f er coefficients. The tubes 
used for data collection were the smooth, modified High 
Flux, and Modified Turbo-3 tubes. The modified Wilson plot 
anaylsis could not be directly run on the enhanced tubes 
since reliable boiling correlations for these surfaces do 
not exist. Also, only two unknowns can be allowed between 
correlations for the inside and outside heat - transfer coef- 
ficients. The steps used in this procedeurs are outlined 
below and are contained in DRPd , • which is cited in Appendix 
A: 

Basic Eouations 

The overall heat - trans fer coefficient in terms of the 
overall thermal resistance is given by: 



o 



(C. 1) 



o 



o 



where 




(C. 2) 



Modified Wilson Plot 
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1. Assume C £ and T f for the boiling side and compute: 

tO . 5 

(C. 3) 





-i j 


/C \ 


t 


o \ 




^Pr 1 ' 7 ) 

Q 


I * | 









„ _2 
U LI i 

~1T 

n c 



( C sf) 



(C.4) 



Note: with r equal to 0.333, equation (C.3) represents the 
Rohsenow correlation [Ref. 23: p. 969]. 

2. Assume C and compute: 



P k - 0 . 3 _ 0 . 33 
r = — Se ?r 




3 . 14 



( C. 5 ) 



h . = c . r 

i i 



(C. 6) 



Substitute hn and h i from steps 1 and 2 above into 
equation (C. I) and rearrange to yield: 

— - R 2 = — 21- + (c (C. 7) 

u ^ c . a . r \ s ' 



4. Now, let: 



(C.3) 



and 



Q.A 



X = 



i c\ . 



(C. 9) 



Construct the least - squares line for Y versus X in the 
f o rm o f : 



Y = mX * C 



( C. 10) 
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5. 



Compute a new set of values 



for C . and C - 
1 3 - 



as follows: 



C. = -i- (C. 11) 

and 

C f = c r (c. 12) 



6. Repeat steps 1 through 5 until convergence for Ci and 
C s f between two successive iterations is less than 0. 1 
percent. 
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APPENDIX D 





DATA RUNS 






The following table outlines the 


data runs used in this 


investigation: 






Where: 




1. 


WH - Smooth tube 




2. 


HF- High Flux ( i. e. , porous - coated Korodense) tube 


3. 


HFM - Modified High Flux (i.e. , 
off) tube 


porous coating machined 


4. 


TB - Turbo-B tube 




5. 


TBM - Modified Turbo-B (i.e. 
machined off) tube 


, external enhancement 


6 . 


T - inlet water temperature 




7. 


V - water velocity- 




8. 


const - constant 




9. 


dec - decreasing 




10. 


inc - increasing 
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TABLE 4 

DATA RUN DESCRIPTION 



File Name Oil. % # points 



Description 



WHO 1 


0 


10 


WHO 2 


0 


10 


WHO 3 


0 


9 


WHO 5 


0 


10 


WHO 6 


0 


8 


WHO 7 


1 


10 


WH08 


1 


10 


WHO 9 


2 


10 


WH10 


2 


10 


WH10A 


2 


10 


WH11 


6 


10 


WH12 


6 


8 


WH13 


6 


8 


HF12 


0 


10 


HF13 


0 


12 


HF14 


0 


12 


HF15 


0 


10 


HF16 


0 


12 


HF22 


0 


7 


HF23 


0 


6 


HF24 


1 


12 


HF25 


1 


10 


HF26 


1 


6 


HF27 


1 


7 


HF28 


2 


12 


HF29 


2 


12 


HF30 


2 


6 


HF31 


2 


8 


HF32 


6 


10 


HF33 


6 


10 


HF34 


6 


7 


HF35 


6 


7 


HFM01 


0 


8 


HFM02 


0 


8 


HFM03 


0 


10 


HFM04 


0 


8 


HFM05 


0 


8 


HFM06 


0 


8 


HFM07 


0 


10 


HFM08 


0 


8 


TBO 1 


0 


10 


TB02 


0 


10 


TB03 


0 


10 



Light Off Effect 
( cold water) 

Light Off Effect 
(warm water) 

Light Off Effect 
earn) 

const V dec 

const V dec 

const V dec 

const V dec 

const V dec 

const V dec 

Repeatability WH10 
T const V dec 

T const V dec 

Light Off Effect 
(warm water) 

T const V dec 

T const V dec 

Light Off Effect 
steam) 

ight Off Effect 
cold water) 
ight Off Effect 
(warm water) 
Repeatability of HF12 
inc V const 
T dec V const 
T const V dec 

T const V dec 

dec V const 

T inc V const 

T const V dec 

T const V dec 

T dec V const 

(Repeatability of HF28) 
inc V const 

T const V dec 

T const V dec 

T dec V const 

T inc V const 

Sieder-Tate Coeff Run 
Sieder-Tate Coeff Run 
Sieder-Tate Coeff Run 
Sieder-Tate Coeff Run 
Sieder-Tate Coeff Run 
Sieder-Tate Coeff Run 
Sieder-Tate Coeff Run 
Sieder-Tate Coeff Run 
Light Off Effect 
( cold water ) 

Light Off Effect 
(warm water) 

Light Off Effect 
( steam) 
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TB04 


0 


10 


TB05 


0 


6 


TB06 


0 


10 


TB07 


0 


7 


TB08 


1 


10 


TB09 


1 


6 


TB10 


1 


10 


TBll 


1 


6 


TB12 


2 


10 


TB13 


2 


10 


TB14 


2 


7 


TB15 


2 


10 


TB16 


2 


6 


TB17 


6 


10 


TB18 


6 


7 


TB19 


6 


7 


TB20 


6 


8 


TB21 


6 


10 


TB22 


6 


10 


TBM04 


0 


8 


TBM05 


0 


8 


TBM06 


0 


6 


TBM07 


0 


6 


TBM08 


0 


8 


TBM09 


0 


8 


TBM10 


0 


8 



T const V dec 

(Repeatability of TB03) 
T dec V const 

T const V dec 

T inc V const 

const V dec 
T dec V const 

T const V dec 

T inc V const 

T const V dec 

T const V dec 

(Repeatability of TB12) 
T dec V const 

T const V dec 

T inc V const 

T const V dec 

T dec V const 

T inc V const 

T const V dec 

ight Off Effect 

steam) 

ight Off Effect 
( cold water ) 

Light Off Effect 
(warm water) 

Sieder-Tate Coeff Run 
Sieder-Tate Coeff Run 
Sieder-Tate Coeff Run 
Sieder-Tate Coeff Run 
Sieder-Tate Coeff Run 
Sieder-Tate Coeff Run 
Sieder-Tate Coeff Run 
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A?? END 17. S 

UNCERTAINTY ANALY SIS 



The uncertainity for mass flow race, Reynolds number, 
heac flux, LMTD , wall resistance, overall heac - transfer 
coefficient, inside heac- transfer coefficient, and outside 
heac - transfer coefficient were analysed for selected runs of 
the smooth, High Flux and Turbo-3 tubes. The analysis was 
based on the Kline-McClintock [Ref. 32: p. 3] method of 
uncertainity analysis. For example, the following equations 
were used for the uncertainity of the heat flux: 



q = 




( 2 . 1 ) 



and 



Q 




( 2 . 2 ) 



In accordance with 
are given by: 




and 



Kline and McClintock, 




the uncertainities 

0 . S 

( E. 3 ) 



60 

Q 



to . /Ai\ 7 • (— 2a- 

j\ a / V = ?/ ^ T Ln* : 3 ur ) 




0.5 

(2.4) 



The uncertainity in the conduction lasses far the unenhanced 
er.d3 of the boiling tube were considered negligible in 
comparison to the boiling surface uncertainity and as such, 
were disregarded. Table 5 lists the uncertainities for the 
previously mentioned values. 
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TABLE 5 





UNCZ2.7AIN7 


"£ ANALYSIS 


? Z3.CZ.IT 


S 




file 


WHO 5 


WHO 5 


Hr 13 


Hr 13 


7304 


1304 


. q 


31140 


26690 


49540 


33730 


6 3 730 


53670 


(w/m ) 


V 


2. 24 


1. 39 


2. 24 


1. 39 


1. 37 


1. 07 


( m/ s ) 


Sift 

A 


1. 56 


2. 51 


1. 5 6 


2. 51 


1. 43 


2. 51 


£Re 


1. 79 


2. 65 


1. 73 


2. 61 


1. 61 


2. 51 


Re 


5o 

q 


1. 5 9 


4 . • -J *m 


1. 5 9 


2. 52 


1. 46 


2. 52 


5i 'm 




4. 96 


3. 59 


3. 10 


2.45 


2. 29 


1. 56 


urra 

53 


* 

R 


2. 76 


2. 76 


11. 17 


11. 17 


2. 92 


2. 92 


W 

u 


5. 21 


■4. 39 


3.43 


3. 52 


2. 71 


3. 02 


3 

5h . 
h. 


2. 03 


2. 56 


1. 61 


2. 25 


1. 43 


2. 21 


i 

5h 

3 

h 


5.34 


5. 03 


4. 90 


4. 39 


3. 35 


3. 96 


3 

Note: 


Ail runs were 


oerf omed 


in pure 


ref riger 


ant 




*rt - 

Hr = 


iraoocn cuDe 
Koroaer.se cube wit 


h High F 


lux Coaci 


ng 




IB = 


-uroa - 3 


tuae 






The la 


rge uncertain! 


cy for Che wall ri 


esis tance 


of f i 



HF13 (porou3- coated Koroder.se Cube) is due Co Che 10 
percar.c uncertainity of Che chemal conductivity coeffi- 
cient. For this investigation, Che coefficient used was an 
average value for Che values given in Che open literature 
for coooer-nickel olates or Cubes. 
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APPENDIX F 

LIST OF NOMENCLATURE 



1. NOMENCLATURE 

A - Surface Area 

c_ - Specific heat at constant pressure 

tP. 

- Oil concentration 

C.; - Inside Sieder-Tate- type coefficient 

- Rohsenow coefficient 

b A 

d^ - Inside diameter of tube 

d Q - Diameter of mouth of reentrant cavity 

D - Diameter 

e - Depth of internal ridge or fin 

E v - Energy transfer per unit volume 

f - Friction factor 

g - Acceleration due to gravity 

g_ - Gravitational constant 

h - Convective heat- transfer coefficient 

hr^ - Specific enthalpy of vapor-liquid mixture 

k - Conductive heat- transfer coefficient 

k m - Conductive heat - transfer coefficient of metal 

L - Tube length 

L e - Entrance effect length 

m - Mass flow rate 

p - Distance between ridge or fin peaks 

Pr - Prandtl number 

q - Heat flux 

Q - Total heat transfer 

r - Tube radiu3 

Re - R.eynolds number 

R.^ - Thermal wall resistance 

St - Stanton number 
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U 0 - 

- 

V - 
A - 

o 

a 

o 

a 

r 

U - 

u o ' 
,J r- 

p - 

P^ - 



Overall heat- trans f er coefficient 
Specific volume of vapor 
Velocity 

Change ( i. e. } T^ n - T cu ^ ) 

Surface tension 

Surface tension of oil 

Surface tension of refrigerant 

Viscosity 

Viscosity of oil 

Viscosity of refrigerant 

Density 

Density of oil 
Density of refrigerant 



2. SUBSCRIPTS 



f - 


fluid 


i 


inside 


in - 


inlet 


m 


mixture 


o 


outside 


o 1 - 


oil- liquid 


out - 


outlet 


ri - 


refrigerant - liquid 


s 


smooth tube 


sat - 


saturation 


w 


water 


wo - 


outside wall 
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